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CHAPTER  ONE 


INTRODUCTION 

Parametrically  excited  multi-degree-of-freedom  systems  have  re- 
ceived considerable  attention.  Linear  models  of  such  systems  predict 
regions  of  instability  in  the  parameter  space  (in  this  case  the  para- 
meters are  essentially  the  frequency  and  the  amplitude  of  the  excitation) 
where  the  amplitude  of  the  response  grows  without  bound.  The  inclusion 
of  damping  in  the  model  does  not  change  the  basic  character  of  the 
resul ts. 

Of  course,  amplitudes  of  real  systems  do  not  grow  without  bound; 
the  amplitudes  are  limited  by  nonlinear  effects.  Moreover,  nonlinear 
effects  make  it  possible  for  a bounded  resonant  motion  to  exist  in 
regions  where  the  linear'models  predict  such  motions  are  impossible. 

The  interaction  of  the  phase  and  amplitude,  which  is  a characteristic  of 
nonlinear  systems,  is  responsible  for  this. 

The  system  of  governing  equations  is  a set  of  second-order  non- 
linear ordinary  differential  equations  having  variable  coefficients.  The 
nonlinearity  considered  in  this  work  is  cubic  in  nature  and  considered 
small.  Thus  the  equations  may  be  referred  to  as  weakly  nonlinear.  The 
excitation  is  considered  to  be  harmonic.  With  some  types  of  parametric 
resonance  several  modes  can  be  strongly  excited  by  a single  harmonic 
frequency.  Internal  resonances  can  be  responsible  for  strong  modal 
coupling  and,  as  a consequence,  for  a significant  transfer  of  energy 
from  one  mode  to  another. 
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In  what  follows  a general  method  of  analyzing  parametrically 
excited  nonlinear  multi -degree-of-freedom  systems  will  be  described  In 
detail  and  applied  to  several  physical  examples.  The  parametric  excitation 
is  assumed  to  be  harmonic  with  constant  amplitude  and  frequency.  Both 
internal  and  parametric  resonances  will  be  considered.  The  resonances  to 
be  examined  include: 

a)  parametric 

\ Z l n.u). 

b)  internal 

to  = ) a. id. 

n L ii 

where  the  n,  n.,  and  a^  are  integers,  X represents  the  frequency  associated 
with  the  excitation,  and  the  oj. , i = 1,  2,  3...,  represent  the  natural 
frequencies  of  the  system.  Modal  damping  is  also  included. 

The  response  of  the  system  is  examined  for  the  case  where  the 
natural  frequencies,  , are  distinct  as  well  as  for  the  case  where 
there  is  one  repeated  frequency,  i.e.  the  so-called  flutter  case. 

Regions  describing  the  existence  of  trivial  and  nontrivial  solutions  are 
discussed  along  with  the  stability  of  these  solutions  within  these 
regions.  Before  proceeding  further  a survey  of  the  literature  concerning 
parametrically  excited  nonlinear  systems  is  in  order. 

1.1  Literature  Review 

A survey  of  the  literature  reveals  that  there  Is  an  abundance  of 
material  written  on  parametric,  nonlinear  and  resonance  phenomenon.  A 
totally  comprehensive  review  would  be  prohibitive  here  both  in  time  and 
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space.  The  most  complete  effort  in  this  direction  has  been  accom- 
plished by  Nayfeh  and  Mook  (1979). 

Faraday  (1831)  appears  to  be  the  first  to  have  observed  the  phe- 
nomenon of  parametric  resonance  in  conjunction  with  surface  waves  in 
fluid-filled  cylinders.  He  was  followed  by  Melde  (1859),  Strutt  (1887), 
Stephenson  (1906)  and  Raman  (1912)  who  worked  with  vibrating  strings. 
Stephenson  (1908)  showed  that  periodic  loading  of  columns  can  have  a 
stabilizing  effect  while  Beliaev  (1924)  showed  that  a lateral  motion  can 
occur  even  though  the  axial  loading  of  a column  is  below  the  static 
buckling  load.  Other  early  works  on  columns  include  Andronov  and 
Leontovich  (1927),  Krylov  and  Bogoliubov  (1935)  and  Chelomei  (1939). 
Numerous  books  on  parametric  excitations  include  Bondarenko  (1936), 
McLachlan  (1947,1950),  Den  Hartog  (1947),  Minorsky  (1947,1962),  Stoker 
(1950),  Hayashi,  (1953a,  1964),  Coddington  and  Levinson  (1955),  Malkin 
(1956),  Cunningham  (1958),  Kauderer  (1958),  Bogoliubov  and  Mitropolsky 
(1961),  Bolotin  (19b4),  Andronov,  Vitt,  and  Khaikin  (1966),  Meirovitch 
(1970),  Cesari  (1971),  Nayfeh  (1973)  and  Evan-Iwanowski  (1976). 

Early  studies  of  nonlinear  vibrations  of  bars  were  conducted  by 
Woinowsky-Krieger  (1950)  and  Burgreen  (1951)  and  involved  the  use  of 
elliptic  functions  in  conjunction  with  an  assumed  single  mode  spatial 
function.  Others  who  investigated  the  free  oscillations  of  beams  with 
hinged  ends  using  various  techniques  include  Wagner  (1965),  Srinivasan 
(1965,1967),  Woodall  (1966),  Evensen  (1968),  Rehfield  (1973,  1975)  and 
Lou  and  Sikarskie  (1925).  Multiple  theoretical  and  experimental  in- 
vestigations were  conducted  by  Eisley  (1964b),  Morris  (1965),  Srinivasan 
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(1966a),  Bennett  and  Rinkel  (1972),  Bert  and  Fisher  (1972)  and  Rehfield 
(1974a),  while  multi-mode  forced  responses  were  investigated  by  Eisley 
and  Bennett  (1970),  Bennett  and  Eisley  (1970)  and  Busby  and  Weingarten 
(1972).  Straight  and  buckled  beams  were  investigated  by  Tseng  and 
Dugundji  (1970,  1971)  using  harmonic  balance.  Others  who  analyzed 
buckled  beams  include  Eisley  (1964a),  and  Min  and  Eisley  (1972). 

Atluri  (1973)  examined  beams  with  large  curvatures  and  Nayfeh 
(1973a)  studied  a beam  with  slowly  varying  properties  along  its  length. 
Both  used  the  method  of  multiple  scales.  Nayfeh,  Mook  and  Lobitz  (1974) 
and  Verma  and  Murthy  (1974)  analyzed  nonuniform  beams  while  Raju, 
Venkateswara  and  Kanaka  Raju  (1976)  studied  tapered  beams. 

A comprehensive  treatment  of  the  vibrations  of  linear  flat  plates 
is  given  in  Leissa  (1969).  ‘Nonlinear  vibrations  of%  plates  has  received 
extensive  treatment.  Yamaki  (1961),  Smith,  Malme  and  Gogos  (1961), 
Eisley  (1964b),  Murthy  and  Sherbourne  (1972),  Bayles,  Lowery  and  Boyd 
(1973),  Vendham  (1975a, b),  Crawford  and  Atluri  (1975),  Venkateswara, 

Raju  and  Kanaka  Raju  (1976a)  all  analyzed  nonlinear  vibrations  of  rec- 
tangular plates.  Anisotropic  rectangular  plates  were  treated  by 
Sathyamoorthy  and  Pandalai  (1970),  Bennett  (1971)  Chandra  and  Basava 
Raju  (1975a, b)  and  Chandra  (1976).  Plates  of  a variety  of  shapes  were 
investigated  by  numerous  authors:  circular  by  Yamaki  (1961),  Nowinskl 
(1962),  Bulkeley  (1963)  Huang  and  Sandman  (1971),  Huang  (1972a,b,  1973, 
1974)  and  Sridhar,  Mook  and  Nayfeh  (1975,  1978);  annular  by  Sandman  and 
Huang  (1971),  Huang  and  Woo  (1973)  and  Huang,  Woo  and  Walker  (1976); 
tringular  by  Vendhan  and  Dhoopar  (1973),  Vendhan  and  Oas  (1975)  and 
Vendhan,  (1975b);  and  elliptic  by  Lobitz,  Nayfeh  and  Mook  (1977). 
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The  phenomenon  of  flutter  in  plates  and  shells  has  received  wide 
attention.  The  book  by  Dowell  (1975)  discusses  the  problem  very  suc- 
cinctly. Other  work  in  this  area  includes  the  study  of  flat  plates  by 
Kobayashi  (1962a,b),  Bolotin  (1963),  Dugundji  (1966),  Dowell  (1966, 

1967a,  1973),  Morino  (1969),  Vetres  and  Dowell  (1970),  Eastep  and  McIntosh 
(1971),  Kuo,  Morino  and  Dugundji  (1972)  and  Smith  and  Morino  (1976). 
Numerous  authors  have  studied  curved  plates  and  shells.  Interaction  of 
panel  flutter  with  parametric  excitation  was  analyzed  by  Dowell  (1970a), 
Dzygadlo  (1970)  and  Kuo,  Morino  and  Dugundji  (1973). 

1 .2  Assessment  of  Previous  Work 

As  demonstrated  in  the  last  section  there  is  a tremendous  volume  of 
material  which  has  been  written  on  free  and  forced  vibration;  parametric 
excitations;  linear  and  nonlinear  systems;  resonance;  and  single,  two- 
degree  and  multi -degree-of-freedom  systems.  An  attempt  to  comment  on 
the  value  and  limitations  of  each  of  these  works  and  how  each  relates  to 
the  current  discussion  would  be  of  marginal  utility.  Instead,  the 
discussion  will  be  limited  to  these  few  recent  works  which  are  the 
most  relavent  and  which  are  significant  contributions  in  themselves. 

Most  authors  cited  here  have  discussed  the  various  aspects  of  the 
analysis  described  in  the  introduction.  However,  no  one  has  attempted 
to  unify  the  analysis  and  make  it  comprehensive  until  now. 

For  instance,  Tseng  and  Dugundji  (1970,  1971)  discussed  straight 


and  buckled  beams  using  harmonic  balance,  but  their  analysis  did  not 
involve  modal  coupling.  Atlurl  (1973)  studied  nonlinear  vibrations  of 
hinged  beams  for  large  curvatures  but  did  not  consider  axial  excitation. 
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Nayfeh  (1973a)  considered  a beam  with  slowly  varying  properties  along 
its  length  but  did  not  consider  modal  coupling.  Yamomoto  and  Saito 
(1970)  analyzed  parametrically  excited  multi -degree-of-freedom  systems 
using  the  method  of  averaging.  However,  the  system  considered  was  linear 
and  did  not  include  repeated  frequencies.  Hsu  (1963,  1965)  also  studied 
parametric  excitation  in  multi -degree-of-freedom  systems  and  developed 
stability  criteria.  He  used  the  method  of  averaging  and  considered 
distinct  frequencies.  Tso  and  Asmis  (1974)  studied  multiple  parametric 
resonances  in  nonlinear  systems.  They  discussed  parametric  resonances 
and  combination  resonances  but  confined  their  analysis  to  two-degree- 
of-freedom  systems.  In  addition,  they  did  not  include  an  all  important 
cubic  term  in  the  analysis.  Fu  and  Nemat-Nasser  (1972,1975)  studied  the 
stability  of  multi-degree-of-freedom  systems  and  included  the  effect  of 
repeated  frequencies.  However,  their  analysis  was  confined  to  linear 
systems.  Another  excellent  work  on  parametrically  excited  linear 
systems  having  many  degrees  of  freedom  was  conducted  by  Nayfeh  and  Mook 
(1977).  Additional  studies  were  conducted  by  Sugiyama,  Fujiwara,  and 
Sekiya  (1968)  using  analog  simulation,  and  they  found  combination 
resonances.  Iwatsubo,  Saigo  and  Sugiyama  (1973)  and  Iwatsubo,  Sugiyama 
and  Ogino  (1974)  performed  theoretical  and  experimental  investigations 
on  parametrically  excited  columns  and  also  demonstrated  the  existence  of 
combination  resonances.  They  included  the  effects  of  internal  and 
external  damping  in  their  analysis  but  not  internal  resonance.  The 
effects  of  Internal  resonance  have  been  virtually  ignored  in  most  works 
until  recently.  The  importance  of  internal  resonance  in  the  analysis  of 
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ship  motions  has  been  demonstrated  by  Nayfeh,  Mook  and  Marshall  (1973, 
1974)  and  Mook,  Marshall  and  Nayfeh  (1974).  Its  importance  in  the 
analysis  of  structural  vibrations  has  been  shown  by  Nayfeh,  Mook  and 
Sridhar  (1973),  Nayfeh,  Mook  and  Lobitz  (1974),  Sridhar,  Nayfeh  and  Mook 
(1975)  and  Mook,  Sridhar  and  Nayfeh  (1978).  These  latter  works  have  been 
a major  advancement  toward  unifying  nonlinear,  resonance  studies. 

However,  they  do  not  include  parametric  excitations. 

The  majority  of  work  discussed  so  far  deals  with  distinct  fre- 
quencies. When  the  natural  frequencies  are  no  longer  distinct  the  response 
enters  the  so-called  flutter  mode  or  condition.  The  book  by  Dowell 
(1975)  provides  an  insight  into  the  mechanism  that  can  cause  flutter. 

Mori  no  (1969)  and  Smith  and  Mori  no  (1976)  have  presented  an  excellent 
discussion  demonstrating  limit  cycles  and  stability.  However,  neither 
included  the  effect  of  an  in-plane  excitation.  The  role  of  a parametric 
excitation  has  thus  been  neglected.  The  numerical  works  of  Dowell 
(1966,  1970)  suffer  from  the  same  deficiency  as  those  of  Kuo,  Morino  and 
Dugundji  (1972).  Constant  iri-plane  loads  were  used  in  these  analyses. 

A harmonic  excitation  and  its  influence  on  flutter  were  considered  by 
Dowell  (1974),  Dzygadlo  (1970)  and  Kuo,  Morino  and  Dugundji  (1973). 
However,  the  excitation  was  not  parametric  in  nature. 

1 .3  Contributions  of  the  Current  Work 

Though  a vast  amount  of  work  has  been  cited  here,  there  are  some 
significant  facets  of  analysis  of  parametrically  excited  multi-degree- 
of- freedom  systems  that  have  not  been  examined.  The  role  that  Internal 
resonance  plays  has  been  virtually  ignored.  Parametric  excitations 
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in  the  nonlinear  flutter  problem  has  been  similarly  neglected.  In 
addition,  there  has  been  no  single  method  of  analysis  that  is  general  in 
nature  and  applicable  to  a wide  variety  of  problems.  This  work  investigates 
these  fascinating  aspects  of  multi-degree-of-freedom  systems  and  thus 
fills  a longstanding  void.  In  the  process  a unified  versatile  approach 
is  developed  and  some  interesting  behavior  is  discussed.  The  following 
contributions  are  unique  and  significant: 

a.  A general  method  for  analyzing  parametrically  excited  non- 
linear multi-degree-of-freedom  systems  is  developed.  To  demonstrate  its 
versatility,  we  apply  the  method  to  the  following  physical  systems: 

1)  Systems  with  distinct  natural  frequencies 

a)  .without  internal  resonance 

b)  with  internal  resonance 

2)  Systems  with  repeated  natural  frequencies 

a)  without  internal  resonance 

b)  with  internal  resonance 

b.  Regions  where  trivial  and  nontrivial  solutions  exist  are 
defined  and  the  stability  of  the  solutions  within  each  region  is  dis- 
cussed. Nontrivial,  unstable  solutions  have  been  shown  to  exist  In 
regions  where  nontrivial  stable  solutions  are  known.  Numerical  solutions 
do  not  hint  at  the  existence  of  these  solutions. 

c.  The  role  of  Internal  resonance  in  parametrically  excited, 
nonlinear  systems  is  explored.  Strong  modal  interaction  is  demonstrated 
as  a consequence  of  the  presence  of  the  cubic  nonlinearity  and  the 
Internal  resonance.  Because  of  this  modal  coupling,  modes  other  than 
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the  one  excited  can  dominate  the  response.  A multiplicity  of  jumps  is 
shown  to  exist. 

d.  Parametric  excitation  in  the  nonlinear  flutter  problem  has 
been  examined  in  detail  for  the  first  time,  and  includes  the  effects  of 
internal  resonance. 

e.  Limit-cycle  behavior  in  the  flutter  problem  is  developed  in  a 
unique  analytic  way.  The  condition  which  predicts  the  onset  of  flutter 
is  developed  by  using  a linear  analysis.  The  solution  grows  without 
bound.  Interestingly,  the  same  condition  in  the  nonlinear  analysis 
predicts  the  existence  of  a real  nontrivial  solution  with  a finite 
amplitude,  the  so-called  limit  cycle. 

f.  The  current  work  is  divided  into  several  parts.  In  Part  I, 
the  analysis  deals  with  distinct  natural  frequencies.  Chapter  2 is 
devoted  to  the  problem  formulation  and  generation  of  the  solvability 
conditions.  Chapter  3 deals  with  the  analysis  of  parametric  resonance 
in  the  absense  of  internal  resonance.  Chapter  4 is  the  heart  of  the 
analysis  and  includes  the  internal  as  well  as  the  parametric  resonance. 

In  Part  II,  the  same  general  approach  is  followed  for  the  analysis 
of  non-distinct  natural  frequencies.  Chapter  5 develops  the  limit-cycle 
behavior  and  the  general  solvability  conditions  for  the  situation  In- 
volving repeated  frequencies.  Chapter  6 considers  parametric  resonance 
in  the  absence  of  internal  resonance  and  Chapter  7 discusses  internal 
and  parametric  resonance. 

Part  III  is  devoted  to  numerical  examples.  Chapter  8 deals  with 
the  case  of  distinct  frequencies  and  uses  the  beam  as  an  example. 

Chapter  9 develops  the  flutter  problem  for  a simply  supported  plate. 


CHAPTER  TWO 


PROBLEM  FORMULATION  AND  METHOD  OF  SOLUTION 

| 

In  this  chapter  the  nonlinear  system  studied  is  defined.  Modifica- 
tion of  the  basic  system  is  made  in  a subsequent  chapter  when  the  case 
of  non-distinct  natural  frequencies  is  discussed.  Solvability  conditions 
are  developed  by  using  the  method  of  multiple  scales. 


2.1  Problem  Formulation 


The  system  to  be  analyzed  is  governed  by  the  following  set  of 
ordinary-differential  equations: 

..  00  00  00  00 

U„  + = e[~  2C0SXt  T fU  - 2CII  + £ l l ? II  II  II  1 

n n n nm  m n n _£■.  _m__U_U_U_J 

m=l  t m=l  p= I q= I nmpq  m p qJ 


for  n « 1,2...  (2.1) 

where  the  un  are  the  distinct  natural  frequencies  corresponding  to  the 

linear  free-osci nation  modes;  e is  a small  dimensionless  parameter;  the 

f are  the  amplitudes  associated  with  the  harmonic  parametric  excita- 
nm 

tion;  X is  the  constant  frequency  of  excitation;  the  Cn  are  the  modal 
viscous  damping  coefficients;  and  the  rnmpq  are  constant  coefficients. 

Experience  has  shown  that  a straightforward  perturbation  expansion 
for  Un,  for  small  e,  results  in  the  emergence  of  secular  terms,  which 
make  this  expansion  not  uniformly  valid  for  large  t.  Hence,  in  order  to 
obtain  a uniformly  valid  expansion,  we  employ  the  method  of  multiple 
scales.  Specifically,  the  derivative-expansion  version  of  the  method  of 
multiple  scales  is  utilized  to  study  the  steady-state  solutions  and 
their  stability. 
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n 


2.2  Method  of  Solution 


Following  the  method  of  multiple  scales  (Nayfeh,  1973),  we  assume 


expansions  for  the  Un  in  the  form 


Un(t;e)  = Uno(To,Ti)  + eUni(T«,Ti)  + . 


where 


To  = t and  Ti  = et 
The  derivatives  become 


(2.2) 


(2.3a) 


= Do  + eDx  + 0(e2) 


(2.3b) 


|p-  = D§  + e2D0Di  + 0(e2) 


where 


Do  = and  Dx  - jyj- 


(2.3c) 


(2.3d) 


Substituting  equations  (2.2)  and  (2.3)  into  equation  (2.1)  and 
equating  the  coefficients  of  e°  and  e,  we  obtain 


Douno  + u)2Uno  - 0 


(2.4) 


D2U  i + ui2U  i = - 2D0D1Uno  - 2cosXT0  l fnmUm  - 2CnD0Un0 
n nn  n nm  m n n 

OO  00  00 

+ l l l rnm  UoUnoUrto  for  n-  1,2,...  (2.5) 

m>l  p«l  q-l  m P q 

We  can  write  the  solution  of  equations  (2.4)  in  the  form 


Uno  3 An(Ti)exp(iwnT0)  + cc 


(2.6) 


for  n * 1,2,..,  where  cc  stands  for  the  complex  conjugate  of  the  pre- 
ceding term  and  the  An  are,  at  this  point,  unspecified  complex  function 


A 
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of  Ti.  They  are  determined  by  eliminating  "troublesome"  terms  from  the 


Substituting  equation  (2.6)  into  equation  (2.5)  leads  to 


D§U  i + (jj2ll  i = - 2ico  (DiA  + C A )exp(ico  To)  - 7 f A {exp("i(co 

0 n1  n n1  nv  1 n n n;  n 0 ' L,  nm  m HL  ' n 

m=  I 


OQ  00  oo 


+ X)T0]  + exp[i(u)_  - X)To]}  + l l 7 r (AAA 

m m=l  p=l  q=l  nmpq  m p q 

exp[i(a)m  + wp  + u>q)To]  + AmApAqexp[i(aim  + <op  - wq)T0] 


+ AmVqexp[i(aJm  ’ “p  + VTo]  + Wqexp[iK  " wp  -VTo]} 


(2.7) 


where  the  overbars  denote  the  complex  conjugates. 


The  terms  containing  the  factor  exp(±ioonT0)  in  equation  (2.7)  lead 
to  so-called  secular  terms  in  Uni,  while  terms  containing  factors  such 
as  exp(inTo),  where  -oon  = 0(e),  lead  to  so-called  small-divisor  terms 
in  U x.  The  former  renders  the  expansion  nonuniform  (Uni/Uno  is  not 
bounded)  as  t increases,  while  the  latter  are  inconsistent  with  the 
assumed  expansion  (2.2)  because  they  raise  eUni  to  the  same  order  as 
Uno.  Consequently,  both  kinds  of  terms  are  "troublesome"  and  must  be 
eliminated  for  a uniform  expansion. 

Secular  terms,  for  example,  are  associated  with  the  combination 

til  + u)_  - co- 
rn P q 

when  m = n and  p * q.  In  addition  to  all  these  possible  combinations, 
there  may  be  small-divisor  terms  associated  with  combinations  having  the 


un  “ awm  + S + cuq  + ean 


(2.8a) 


where  a,  b,  and  c are  integers  such  that  |a|  + |b|  + |c|  = 3,  and  is 
a detuning  parameter.  If  equation  (2.8a)  is  satisfied,  an  internal 
resonance  is  said  to  exist.  Finally,  if  there  are  combinations  having 


the  form 


X = a ui  + bw  + cp 

m n 


(2.8b) 


where  a and  b are  integers  such  that  |a|  + |b|  = 2,  and  p is  another 
detuning  parameter,  additional  small-divisor  terms  arise.  If  equation 
(2.8b)  is  satisfied,  a parametric  resonance  is  said  to  exist. 

It  follows  that  the  general  form  of  the  conditions  that  eliminate 
the  troublesome  terms  from  the  Uni  is 


21“n<°'\  * CnV  + An  pI,  Vp\  * ‘n  * En  * 0 for  " = ’•2— 

(2.91 


where 


if  p = n 


2^rnnpp  + rnpnp  + rnppn^  if  p ^ n 


The  terms  In(A!,A2,  ...)  and  En (A  x ,A2 , . . . ) are  complex  functions,  re- 
sulting from  the  internal  and  parametric  resonances,  respectively,  if 
any  exist. 


2.3  Summarv 


A general  approach  to  the  solution  of  equations  (2.1)  has  been 
developed  using  the  method  of  multiple  scales  to  obtain  a uniformly 
valid  solution.  The  solvability  conditions  yielded  equations  (2.9), 
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which  are  general  and  can  be  applied  to  a large  variety  of  problems, 
depending  on  the  form  of  the  internal  and  parametric  resonances  in- 
volved. 


CHAPTER  THREE 

PARAMETRIC  RESONANCE  IN  THE  ABSENCE  OF  INTERNAL  RESONANCE 

In  this  chapter  parametric  resonance  is  examined  in  detail.  The 
effect  of  an  internal  resonance  is  discussed  in  the  next  chapter. 

Various  cases  are  considered  and  some  general  observations  are  given. 

The  emphasis  is  on  the  steady-state  response  and  its  stability. 

3. 1 No  Internal  Resonance 

In  this  case  I s 0 for  all  modes  and  E„  = 0 for  all  modes  that  are 
n n 

not  part  of  any  resonance  with  the  excitation.  The  Ep  for  the  resonating 
modes  depend  on  the  resonant  combinations  of  frequencies.  In  the  first 
case  considered,  there  is  neither  internal  nor  parametric  resonance. 

3.2  The  Case  of  Modes  Not  Involved  in  a Parametric  Resonance 

For  modes  not  associated  with  a parametric  or  internal  resonance 


I = E i 0 in  (2.9).  Now  we  let 
n n 


An  ■ \ an(Ti)exp[ian(Ti)] 


(3.1) 


where  the  amplitude  an  and  the  phase  an  are  real  functions  of  Ti.  Upon 
substituting  (3.1)  into  the  solvability  condition  (2.9)  and  separating 
the  resulting  equation  into  its  real  and  imaginary  parts,  we  obtain 


u>_(a'  + C a ) = 0 
n n n n 


+lanl  Wn  - 0 

P 


'n_n”n  8 n t 'np~p 


(3.2a) 

(3.2b) 
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where  the  primes  designate  differentiation  with  respect  to  Tj.  The 
solution  of  (3.2a)  is 


an  = a0exp[-CnTi] 


(3.3) 


and  it  follows  that  an  + 0 as  Ti  -*■  »;  so  that  the  steady-state  amplitude 


an  = 0.  Therefore,  the  steady-state  solution  has  the  form 
un  ■ °<£> 


(3.4) 


Consequently,  only  the  directly  excited  modes  and  the  modes  involved  in 
an  internal  resonance  can  be  part  of  the  first  approximation  of  the 
steady-state  response. 

3.3  The  Case  Where  X T 2wk,  n f k 

The  first  case  to  be  considered  is  the  case  of  a parametric  reson- 
ance where  X ~ 2u)k.  A detuning  parameter,  p,  is  introduced  and  used  to 
express  the  nearness  of  X to  2wk  quantitatively  according  to 

X 3 2o>k  + ep  (3.5a) 

It  follows  that 


(X  - u)k)T0  3 oj^To  + eToP  * oj^Tq  + pTi 


Then  for  n f k the  resonant  terms  become 


En  3 0 for  n + k 


(3.5b) 


(3.5c) 


and 


Ek  * * ^ kk^kexP( i oT 1 ) * — a|cexPC‘* (pT i - ak)]  (3.5d) 
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Substituting  (3.5)  into  the  solvability  conditions  (2.9)  and  separating 
into  real  and  imaginary  parts,  we  obtain 


kk 


and 


uk^ak  + ^kak^  + ~T~  aks^nii  = ® 


1 ® fkk 

Vkak  + ak  i}  Ykpap  ‘ ~T  akC0Su  = 0 


(3.6a) 

(3.6b) 


where 


y = pTi  - 2afc 


(3.6c) 


For  the  steady-state  solution  all  the  3^=0  and  y'  = 0.  From  equations 

(3.3)  and  (3.4)  it  follows  that  an  - 0 as  Ti  -*•  » so  that  the  steady- 

state  amplitude  an  = 0 for  n i1  k.  The  trivial  solution  ak  = 0 is  a 

solution  of  equations  (3.6).  For  nontrivial  solutions  ak  f 0,  equations 

% 

(3.6)  reduce  to 

wk^k  + — siny  - 0 (3.7a) 

and 

"F  + S'  Ykkak  ‘ ~7T  COSu  = 0 (3-7b) 

where 

y1  = 0 = p - 2ak  (3.7c) 


Equations  (3.7)  may  now  be  solved  for  ak  and  y. 
(3.7a)  and  (3.7b)  leads  to 

ir2  + + 1 v a2!2  - 

ukck  + + S YkkV 


Squaring  and  adding 


(3. 7d) 
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We  note  that  for  nontrivial  solutions  to  exist 


lfkk'  - 2ukCk 

when  p > 0 and  ykk  < 0.  When  p < 0 and  ykk  < 0 


f2 


P2  < — r-  - 4C  2 


OJ 


(3. 7e) 


(3.7f ) 


is  required  for  nontrivial  solutions.  A more  complete  discussion  of  the 
regions  where  solutions  exist  is  contained  in  the  numerical  examples  in 
Chapter  8.  We  note  here  that  ykk  < 0 in  those  examples. 

When  equation  (3.7d)  is  substituted  into  equations  (2.2)  and  (2.6), 
the  resulting  steady-state  solution  becomes 


Un  = 0(e)  n * k 


and 


(3. 7g) 

(3. 7h) 


Uk  = akcos[(cuk  + e £)t  + tk]  + 0(e) 

where  the  phase  xk  is  a constant  that  depends  on  the  initial  conditions. 

The  stability  of  the  steady-state  solution  can  now  be  studied  by 
determining  the  behavior  of  the  system  when  it  is  perturbed  slightly. 
Thus,  we  let 


ak  * ak  + Aak 


U + Au 


(3.8a) 

(3.8b) 


where  ak  and  u represent  steady-state  values.  Substituting  (3.8)  into 
(3.6)  yields 

2wk(a£  + Aak  + Ckak  + CkAak)  + fkk(ak  + Aak)sin(y  + Ay)  * 0 

(3.9a) 


uik(ak  + Aak)otk  + g-  YkkUk  + Aak) 3 - -jp  (ak  + Aak)cos(y  + Ay)  = 0 

(3.9b) 

After  expanding,  neglecting  products  of  small  terms,  and  noting  that 

A 

ak  = 0,  we  obtain 

A A A A 

2o>k(Aak  + CkAak)  + 2u)kakCk  + fk|(Aaksiny  + fk|<aksinu 

A A 

+ fkkakcosyAy  =0  (3.9c) 


ukak  {QZ^L~)  + ujkAak(f)  + s Ykkak + 1 YkkakAak " Hr  akcosy 


fkk  ~ * f kk  /' 

+ ^ 1 aksinyAy  - — Aakcosy  = 0 

Recalling  the  steady-state  equations  and  letting 

(Aak,Ay)  « exp(nTi) 

we  obtain 

A A 

(2u>kft)Aak  + (fkkakcosy)Ay  * 0 

' (T  Ykkak)Aak  + (Hr  + u,kCk)Au  = 0 
from  which 

“kfl!  + 2“kCkn  * T *kk»kfkkcos“  ' 0 


fl2  + 2Ckn  + l^jr"  cosy  * 0 


(3.9d) 


(3.10a) 


(3.10b) 


(3.10c) 


(3.1 Od ) 


which  has  the  solution 


a - - C.,  ±V  C?  - 


f,, 


k - “TG* 


(3.11) 
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For  stable  solutions  the  real  part  of  n must  be  negative,  or  zero  i.e.. 


Re(fl)  < 0 (3.12) 

or  for  stability 

A 

cos  u _>  0 (3.13) 

Hence,  the  p for  stability  can  be  determined  from  equation  (3.7d). 
Additional  discussion  of  the  ramifications  of  the  stability  conditions 
will  be  presented  in  Chapter  8 in  the  numerical  examples. 

3.4  The  Case  Where  X - com  + m + k 

The  case  of  parametric  resonance  where  X ~ + uj^  is  analyzed. 

The  detuning  parameter  p is  now  used  to  measure  the  nearness  of  X to  ^ 
according  to: 


X = + + ep 

(3.14a) 

It  follows  that 

(X  - ) To  3 u^To  + er0p  3 ^To 

+ pT  i 

(3.14b) 

and 

(X  - aj^)To  3 <y0  + eT0p  3 uimTo 

+ pTi 

(3.14c) 

E 3 0 for  n / k or  m 
n 

(3.15a) 

Em«Pd^T.)  - - f^exptMx-^lT.]  . . f^expCK^T,  ♦ pT,)] 

(3.15b) 

and 

Ekexp((»kT.)  ■ - » - fkmV*P[1<"kT*  * °T'>J 

(3.15c) 


i 


21 


Using  the  polar  notation  in  equation  (3.1),  we  can  write 

Em  = ’ “It1  expCi(pTl-ctk)] 
and 


(3. 15d) 


Substituting  equations  (3.15)  into  equation  (2.9),  dividing  by  exp(ictn), 
and  separating  the  resulting  equations  into  real  and  imaginary  parts,  we 
obtain 


^—3  ) + ~ 7“  9|.Siny  = 0 
m'  m mm  2 k 


wk(ak  + Ckak)  + ~1T  amsiny  = 0 

oo  ^ 

iamam  + S’  am  Ympap  * TT  akCOSu  = 0 

Vk°k + 1 ak  1 Vp  - "T  amC0Sy  s 0 


oj. 


where 


(3.16a) 

(3.16b) 

(3.16c) 

(3.1 6d ) 


y = pTl  " ak  - am 


(3. 16e) 


For  n f k,  or  m,  the  amplitudes  and  phases  are  governed  by  equations 
(3.2)  and  (3.3);  thus,  these  modes  do  not  appear  in  the  first  approxima- 
tion. For  the  steady-state  solution  all  = 0 and  u‘  = 0.  It  follows 
that  a^  * am  = 0 is  a solution  to  equation  (3.16),  but  nontrivial  solution 
(i.e.,  ak  t 0 and  am  f 0)  are  desired.  Hence,  the  steady-state  version 
of  equations  (3.16a)  and  (3.16b)  becomes 


wmCmam  + ~T  aksiny  = 0 
P , + fkm  a . a Q 

“kCkak  + ~r  amsinw  0 


(3.17a) 

(3.17b) 


Solving  equations  (3.16c,d)  for  and  a£,  respectively,  and  substituting 
the  results  into  equation  (3.16e),  when  differentiated,  leads  to 
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m 


m 


\ 

(3.18a) 


i_«m  a l^l, 

mm  k k 


Combining  equations  (3.17a)  and  (3.17b)  leads  to 


m 


_ ,/  “kVmk 
~v 


m m km 


(3.18b) 


Hence,  for  nontrivial  solutions  to  exist  and  fkm  must  have  the  same 
sign.  Substituting  equation  (3.18b)  into  equation  (3.18a)  and  dividing 
by  a£,  we  obtain 

»<  Sfe )1/2  - ? sr  (1  * h™  * ff  < Sot  >‘/! « 


m m km 


[(  + ^m)(  ) + ( JK.  + -HL)]aj|  = 0 (3.18c) 

m k m m km  m k 


Yi 


kmx 


"kCkfmk 


Ymk  . Ykkx 


and  hence 


ak  * 


/>  Is*  (1  * 

a/  7“m 


Ck  x__  UJkCkfmk  x 7z 

ir  Kosv  “ p(  rr7T } 

m m m km 


/l  /‘V'k^mk  xl/2  r/  Ymm ^ Ykm  x^k^k^mk  x . Ymk  , Ykk  xi 

c “T  “T  nT  “T  1 


(3.1 8d ) 


The  steady  state  solution  now  has  the  following  form: 


Un  * 0(e) 


n f k or  m (3.19a) 

Um  * amcosC(a)m  + “hi^  + Tm]  + 0(e)  (3J9b) 

* aj(cos[(c0|(  + ea£)t  + t^]  + 0(e)  (3.19c) 

where  tm  and  are  constants  depending  on  Initial  conditions.  The 

nonlinearity  adjusts  the  frequencies  so  that 
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“m  + eot^  + uk  + eak  = % + wk  + £p  = X (3’20) 

If  both  solutions  (am  = ak  ■ 0 and  am  ^ 0,  ak  t 0)  are  stable,  then  the 
initial  conditions  determine  which  solution  represents  the  response. 

The  stability  of  the  system  may  be  examined  by  the  same  method  used  in 
the  case  of  X near  2uk- 

3.5  The  Case  Where  X Z wk  - um>  m ^ k 

The  case  of  parametric  resonance  where  X ~ ajk  -odm  is  analyzed.  The 
detuning  parameter,  p,  is  introduced  and  used  to  express  the  nearness  of 


X to  - ojm , according  to 

X = a)k  - o^m  + ep 

• from  which  we  can  write 

Uk  - X)T0  = (u^  - ep)T0  = 


t±)mT  o ■ pT  i 


(X  + wm)To  = (wk  + ep)T0  = cJkT0  + pTi 


En  = 0 for  n / k or  m 


(3.21a) 


(3.21b) 


(3.21c) 


(3.22a) 


EkexpC1wkT„]  ■ - fkn,A,„exptHx  * ^)T,]  ■ - fkn,A„,«p[1 <“kT"  * » T‘» 

(3.22b) 

EmexpCI^T,]  ■ - fmkAk«p[1(%  - X)T,]  • - VVxpHI^T,  - pT,)] 

(3.22c) 


Using  the  polar  notation  in  equation  (3.1),  we  can  write 


^kmam 


exp[i (pTi  + c^)] 


(3. 22d) 
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and 


Em 


-^jr~  exp[i  C-pTi  + ak)] 


(3.22e) 


Substituting  equation  (3.22)  into  equation  (2.9),  dividing  by 
exp(iak)  and  exp(iam),  respectively,  and  separating  the  resulting 
equations  into  real  and  imaginary  parts  leads  to 


where 


%lam  + 

Snam^ 

mk  k 
~7~ 

siny  = 0 

(3.23a) 

“k(ak  + 

Ckak>  + 

^kmam 

~~r~ 

siny  = 0 

(3.23b) 

<i)mamotm 

a 00 

Ympap 

fmkak  _ n 

~2 — cosy  = 0 

(3.23c) 

“k*kak 

l 

av  00 

Ykpap 

f,  a 

Km  m — a 

2 — COSy  = 0 

(3.23d) 

y * pTr 

-ak  + 

am 

(3.23e) 

For  n t k or  m,  the  amplitudes  and  phases  are  governed  by  equations 
(3.2)  and  (3.3).  For  the  steady-state  solution  all  a^  ■ 0 and  y'  = 0. 


It  follows  that  ak  * am  = 0 is  a solution  to  equations  (3.23),  but  non- 


trivial solutions  (i.e.,  ak  t 0 and  am  / 0)  are  desired.  Hence,  the 
steady-state  version  of  equations  (3.23)  becomes 


uCa  - siny  » 0 

m m m c 


“kck»k + s1nu  ,0 


(3.24a) 

(3.24b) 


and 


1 + Ikm-ja2  + 1 j-£mk  + M]a2 

8 ufo  u>k  m o *■'  K 


“k  ■ k 


1 ♦ i^]cosu  • 0 


nrm 


ak“k 


(3.24c) 


f 
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which  leads  to 


^Cbf, 


k k'mk 


(3.25) 


Thus,  a solution  can  exist  only  If  fkm  and  f^  have  the  different  signs. 
In  this  case,  the  result  for  ak  can  be  obtained  from  (3.1 8d ) by  simply 


changing  the  sign  of  oJm. 


3.6  Summary 


Various  cases  concerning  a parametric  resonance  have  been  discussed. 
The  following  observations  are  pertinent  in  the  absence  of  an  Internal 
resonance: 

(a)  Only  modes  which  are  directly  excited  by  a parametric  resonance 
are  part  of  the  steady-state  response. 

(b)  For  combi  nation- type  resonances  the  lower  mode  appears  to  be 
dominant  in  the  response. 

(c)  The  steady-state  response  is  essentially  that  of  the  forced 
linear-oscillation  modes  with  the  frequency  being  adjusted. 

(d)  Stability  criteria  have  been  outlined  when  more  than  one 
stable  solution  exists,  the  initial  conditions  will  determine 
the  response. 


More  insight  into  the  response  in  the  presence  of  parametric  resonance 
will  be  gained  in  Chapter  8 where  specific  numerical  examples  are 
presented. 


CHAPTER  FOUR 


THE  EFFECT  OF  INTERNAL  RESONANCE 

In  this  chapter  we  discuss  the  effect  of  an  internal  resonance. 
Specifically,  we  examine  in  detail  the  case  where  ~ 3u>m.  Part  of  the 
rationale  for  this  choice  is  that  this  combination  occurs  in  physical 
systems  such  as  the  beam  considered  in  the  numerical  example  in  Chapter 
8. 

The  parametric  resonances  of  the  last  chapter  are  examined  in  the 
presence  of  this  internal  resonance.  See  Chapter  2 for  a review  of 
internal  resonance. 

An  internal  resonance  emerges  when  equation  (2.8a)  is  satisfied. 
When  the  form  of  equation  ?2.8a)  is 

U)k  = 0^  + (dp  + Wq  + eo  (4.1 ) 

it  follows  from  equation  (2.7)  that  the  terms  involved  in  the  internal 


resonance 

are  given  by  the  following: 

h ■ 

h AmApAqeXp(‘ioTl) 

(4.2a) 

V 

rmAkApVXP(i°Tl) 

(4.2b) 

V 

VkVqexp(1oT|) 

(4.2c) 

V 

?qWpeXp(1aTl) 

(4. 2d) 

A 

r = 

k 

rkmpq  + rkpmq  + rkpqm  + rkqpm  + rkqmp  + rkmqp 

(4.3a) 

/\ 

T * 
m 

r + r + r + r 4*r 

mkpq  mkqp  mpkq  mqkp  mpqk  mpqk 

(4.3b) 
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A 

r = r + r + r +r 

p pmqk  ‘pqmk  pmkq  pqkm 

+ rpkmq  + rpkqm^‘3c^ 

A 

rq  ' rqmpk  + ‘rqpmk  + *rqmkp  + rqpkm 

+ rqkmp  + rqkpm^'3c^ 

When  the  form  of  equation  (2.8a)  is 

oj.  = oj„  + 2u_  + ea 
k m p 

(4.4) 

it  follows  from  equation  (2.7)  that 

'k  ’ VmVexpMoTl) 

(4.5a) 

!m  ’ VkVeXp(1oT,) 

(4.5b) 

A 

I = r A.  AA  exp(iaTi) 

p p k m p r ' 

(4.5c) 

/V 

r = r + r + r 

k kmpp  kpmp  kppm 

(4.6a) 

A 

r = r + r + r 

m mppk  mpkp  mkpp 

(4.6b) 

A 

r~r  + r + r +r  + r 

p pmpk  ppmk  ppkm  pmkp 

pkpm  + rpkmp  (4.6c) 

When  the  form  of  equation  (2.8a)  is 

“k  = 3tom  + eo 

(4.7) 

it  follows  that 

‘k  ■ ?k\’«p(-i°T1) 

(4.8a) 

>m  * ?mAkIm2exp(,c,T‘) 

(4.8b) 

where 

rk  ” rkmmm 

(4.9a) 

L 
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r s r + r + r 
m mmmk  mmkm  mkmm 

When  the  form  of  equation  (2.8a)  is 

oj.  = w + i»j  - ui  + ea 
k m p q 

it  follows  that 

!k  ■ wPV‘pHoTi) 

A 

I = r A.  A A exp(iaTi) 
m m k p q 

/\ 

In  = f — A* A_A  exp ( i oT i ) 
p p k m q r 

/\ 

I.  = r A A A.  exp(-ioTi) 
q q m p k r ' ' 

AAA  A 

where  r^,  rm,  rp,  and  Fq  are  given  by  equations  (4.3). 

Other  possible  internal  resonances  can  be  treated  in 

4.1  The  Case  Where  X ~ 2w  and  u>.  ~ 3u  . 
m k m 

The  first  case  to  be  considered  is  X ~ 2w  , w.  ~ 3u>  . 

m k m 

parameters  are  defined  in  the  following  equations. 
wk  = + ea 

X = 2a)m  + ep  = u>k  - ii^  + e(p  - a) 

Em  ■ * fmkVxp[-  Ko  - a)!-,] 

’ ,™Vxp(,°T‘) 


(4.9b) 

(4.10) 

(4.11a) 

(4.11b) 

(4.11c) 

(4. lid) 

a similar  way. 

The  detuning 

(4.12a) 

(4.12b) 

(4.13a) 


Ek  ’ ’ fkAexpCi(c  - °>T'l 


(4.13b) 
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I_  = I’mAtAm2exp(iaTi) 
m m k m r' 


(4.14a) 


!k  = rkAm3exP(‘iaTl)  (4.14b) 

Using  the  polar  form  (3.1),  substituting  (4.13)  and  (4.14)  into  equation 
(2.9),  neglecting  all  except  the  m-th  and  the  k-th  modes,  and  separating 
the  results  into  real  and  imaginary  parts,  we  obtain 

A 


+ W ' / akam2siny  + 


^mmam 


~Sr^  a^sin(6  - u)  = 0 


(4.15a) 


VA'  + I am{Wm2  + Ymkak2)  + fT  akam  C0Sy 
fa  f , 

- JMJJH  cos6  . _mk  ^ cos(g  - y)  = 0 (4.15b) 

/v 

ric 

U)k(ak’  + Ckak^  + “t  am3-sinu  + — sin^e  - u)  = 0 (4.16a) 


“kVk'  + / (l'kmam!  + \kak>  * T am  cos“ 


^kmam 


cos(S  - u)  = 0 


where 


w = aTl  ‘ 3am  + ak 


6 " pTl  ’ 2am 


(4.16b) 


(4.17a) 


(4.17b) 


For  the  steady-state  solution  it  follows  that  a = a,  * 0 is  a solution. 

m k 

For  nontrivial  solutions 


“'■"'V  - 3cSn'  = 0 


(4.18a) 
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*'  = p - 2<V  = 0 


(4.18b) 


which  leads  to 


< ■ f *"d  v ■ ¥ 


Equations  (4.16)  become 


(4.19) 


CO. 


fa 
mm  m 


mk 


1 


mCmam  + ~2  si"6  ' “T~  aksin(B  - v)  - f Vkam2siny  = 0 

(4.20a) 


oa  y 
H m . 'mm  3 
o>  -5-  + — 5—  a"  + 
m2  8 m 


Ymkam  2 'mm~m 


f a. 


8“  ak  - PP  C0S6  - ^ C0S<B  : u) 


* S Ww  cosu  ' 0 


“kCkak  * ^r1  s1n(S  - «>  + 5 fkams1nu  ’ 0 


(4.20b) 

(4.20c) 


fL„a. 


/3o  1 , 'mk  . .2  , Ykk  .3  km  m \ 

ak  <r  - o)  * ~r  Vm  * -r  ak  - ~r- cos(s  • u> 

* 5 fkamcos“  10 


(4.20d) 


A detailed  discussion  of  the  behavior  of  the  solution  of  equations 
(4.20)  is  given  in  Chapter  8 in  the  numerical  examples. 


4.2  The  Case  Where  X - 2u>k  and  u>k  ~ 3a)m. 


The  next  case  to  be  considered  is  X ~ 2u>k,  wk  ~ 3^.  The  detuning 
parameters  are  defined  as  follows: 

X * 2a j.  + ep  (4.22a) 


J 


X 2 2o^  + ep 
o>k  r 30^  + ea 


(4.22b) 
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The  parametric  resonant  terms  are  of  the  form 


E = 0 (4.22c) 

m 

Ek  = - fkkAkexp(ipTi)  (4.22d) 

and  the  internal  resonant  terms  are  the  same  as  those  in  equations 
(4.14).  Substituting  (4.22cd)  and  (4.14)  into  equation  (2.9)  and 
separating  the  result  into  real  and  imaginary  parts,  we  find  that 


%,<a; 

+ C a ) 
m m 

1 A 
4rm 

akams’nu  = 0 

(4.23a) 

(inamam 

a 

^ Ymmam  H 

1 A 

h y , a,2)  + 7r  r a.  a2cosu  3 0 
'mk  k'  8 m k m M 

(4.23b) 

“kK  * 

ckak> 

. fkkak 
+ 2 

sing  + ^ rka^sinu  = 0 

(4.23c) 

Vkak 

* r (Wm  + 

fj-J-Sl.  -I  /v 

ykkak)  - T" cos6  + a rkamcosu 

3 0 

(4.23d) 

where 


6 = pTi  - 2ak  (4.24a) 

y = oTi  + ak  - 3am  (4.24b) 

For  steady  state  solutions,  it  is  easily  seen  that  a = ab  3 0 is  a 

m k 

solution.  For  nontrivial  solutions  equations  (4.24)  become 

S*  = p - 2ctk’  3 0 (4.25a) 

p'  = a + ak'  - 30^'  3 0 (4.25b) 


so  that 

v ■! 


“In'  ’ 


and 


3 


(4.26a) 
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and  equations  (4.23)  become 


oj  C - i r a.  a sinu  = 0 
mm  b m k m 


0 + § . Ymmam  . Ymkak  1 2 _ n 

-T~  + 8 + —5~  + 7 rm  W0Sy  - 0 


'\Ckak  + 


fkkak 


I?  a 3 , 


sine  + g r^a^siny  = 0 


(4.27a) 


(4.27b) 


(4.27c) 


, P . Ykm  , 2 . Ykk  , 

Vk  7 + T akam  + ~T  ak  • 
* 1 Vmcosu  ‘ 0 


fkkak 


(4.27d) 


A detailed  discussion  of  the  behavior  of  the  solutions  of  equations 
(4.27)  is  given  in  Chapter  8 in  the  numerical  examples. 

It  should  bq  noted  here  that  it  is  possible  for  am  - 0 and  a^  f 0, 
in  which  case  equations  (4.27c,d)  become 


ii)  r . 1 kk  k _ n 

kCkak  2 sin  " 0 


U)  _ £ . Ykk  , fkkak  . n 

kak  7 + ~r  ak  • “i-  cos8  0 


(4.28a) 


(4.28b) 


This  is  exactly  the  same  type  of  response  that  occurs  when  there  is  no 
internal  resonance,  equations  (3.7a,b).  Thus,  the  only  mode  which 
responds  is  the  one  which  Is  directly  excited.  This  is  in  contrast  to 
the  case  where  the  frequency  of  the  excitation  X is  near  2uim  and  a 
single-mode  response  is  not  possible. 
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4.3  The  Case  Where  X ~ 4u>m  and  ~ 3o)m. 

The  last  case  to  be  considered  is  for  X ~ 4^  and  ui^  - 3a)m. 

In  the  presence  of  this  internal  resonance  the  parametric  resonance  can 
be  expressed  equally  well  as  X ~ u + co^.  We  may  write 

X = 4u^  + ep  = o)m  + + e(p  - a)  (4.29a) 

where 


OJ.  = 3d)  + ec 
x m 

The  terms  involved  in  the  parametric  resonance  are 


Em  = - fmk\exp[1(o  - a)T'] 


h ’ ' fk,A,expCi(p  -a>T‘] 


(4.29b) 

(4.30a) 

(4.30b) 

(2.9), 
into  real 


Substituting  equations  (4.30)  and  (4.14)  into  equation 
introducing  the  polar  form  (3.1),  and  separating  the  result 
and  imaginary  parts,  we  obtain 

co(a'  + C a ) + m-iv-  sinS  - i T a.a^sinu  * 0 (4.31a) 

mm  mm  Z o m K m 

f ..a, 


Vlflil  * 3 \ + Ymkak>  * C0SS 

* 5 akamc0Sli  ‘ 0 

“k(ak  + Ckak>  + sin8  * 3 fkans1np  * 0 

Vk“k  ♦ r (Yh,*;  * vkkak>  - -t1  cos8 


* J amCOSu  ’ ° 


(4.31b) 

(4.31c) 


(4.31d) 


where 


e = (p  - a)  Ti  - ak  - am 


(4.32a) 


y = aTi  + ak  - 30^  (4.32b) 

For  the  steady-state  solution,  it  follows  that  am  = ak  = 0.  For  the 


nontrivial  solution. 


B'  * » - 0 - V - V = 0 


(4.33a) 


y'  = a + ak'  - 30^'  « 0 


(4.33b) 


which  lead  to 


■ | a"<* 


(4.33c) 


The  steady-state  equations  can  be  written  as 


■Wm  * -E!rL  s1n6  - J "r,Vms(nu  * 0 


(4.34a) 


. p , amYmm  , amakYmk  *mkak  „ 

Vmt  +-8-  + “l 5“  cose 


+ I rm  akamC0Sy  = 0 


(4.34b) 


“kckak  + s("6  * S Vm*1’"1  ' 0 


(4.34c) 


3 2 

3p  - 4a,  akYkk  akamYkm  fkma 

— CT + x + x - X” 


+ s rk  <cos»  = 0 


(4.34d) 


A detailed  discussion  of  the  response  predicted  by  these  equations 


is  given  in  Chapter  8. 


b 


I 
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4.4  Summary 

Internal -resonance  terms  involving  ~ 3wm  have  been  developed. 

The  response  in  the  presence  of  this  internal  resonance  and  several 
parametric  resonances  has  been  examined.  The  following  comments  are 
pertinent: 

a.  It  has  been  shown  that,  even  in  the  presence  of  internal 
resonance,  it  is  possible  for  only  the  mode  which  is  excited  to  enter 
the  response.  The  behavior  is  then  essentially  identical  to  that  in  the 
absence  of  internal  resonance. 

b.  In  all  cases  discussed,  strong  modal  interactions  due  to  the 
presence  of  internal  resonance  are  possible.  A full  discussion  of  the 
dominance  of  one  mode  over  the  others  is  presented  in  Chapter  8,  which 
contains  the  numerical  examples. 


CHAPTER  FIVE 


THE  ANALYSIS  OF  SYSTEMS  HAVING  REPEATED  FREQUENCIES  IN  THE 
ABSENCE  OF  INTERNAL  AND  PARAMETRIC  RESONANCES 

In  the  last  three  chapters  we  have  discussed  the  response  of  a 
nonlinear  system  subject  to  a parametric  excitation.  Throuqhout  the 
entire  analysis  we  have  considered  the  natural  frequencies  to  be  dis- 
tinct. This  is  valid  for  the  beam-column  we  discuss  in  Chapter  8 as 
well  as  for  many  other  physical  systems.  The  analysis  contained  in  the 
next  three  chapters  deals  with  the  situation  where  one  of  the  natural 
frequencies  is  repeated.  As  an  example,  this  occurs  in  plates  subject 
to  aerodynamic  loading  at  the  onset  of  flutter.  Two  natural  structural 

frequencies  merge  due  to  aeroelastic  coupling.  For  a discussion  of  the 

• 

mechanism  of  flutter  the  reader  is  referred  to  Dowell  (1975). 

In  this  chapter  we  follow  a format  similar  to  that  of  Chapter  2. 

The  governing  equations  (2.1)  are  modified  to  account  for  the  repeated 

frequency  and  the  aerodynamic  loading.  In  addition,  scaling  factors 

6i 

will  be  introduced  in  the  form  e to  insure  that  terms  interact  at  the 
proper  order.  The  method  of  multiple  scales  is  utilized  to  generate  the 
solvability  conditions.  The  conditions  that  predict  the  onset  of 
flutter  are  developed  and  limit-cycle  behavior  Is  demonstrated.  We  begin 
now  with  the  governing  equations  when  one  natural  frequency,  u. , is 

J 

repeated,  thus  we  write 
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U_  + u>X  = e F[-(2  cos  At)  l f UJ 
n n n nm  in 

m=  I 

8~  6W  °°  00  «> 

'[2C  0n]  + e \ l l l rnmnU  u u ] 

n n m=l  p=]  q=i  ™pq  m p qJ 


- e 


+ e A[  A I 8UJ  n 3 1,2...  n f k (5.1a) 
1 nm  r 


m=l 


and 


6C  00 
r r 


Uk  * »juk  ♦ Uj  = s r[-(2  MS  Xt)  fkmum: 


5 5 30  00  00 

- e C[2CUUJ  + e N[  l l [ r,  „ UUUl 

k k m*l  p=i  q=i  kmPP  m P P 


5.  » 

+ e [A  l efcmU.  ] where  k * j + 1 


m=l 


(5.1b) 


where  k s j + 1,  A is  a detuning  parameter  (associated  with  the  aero- 
dynamic  loading)  to  be  explained  later  in  a numerical  example,  and 
6p,  <5q,  6n,  and  6^  are  constants  to  be  chosen  so  as  to  insure  the  proper 
interaction  for  the  various  resonances  to  be  considered. 


We  note  that,  if  the  nonlinear,  damping  and  forcing  terms  were 
zero,  the  k-th  mode  would  grow  indefinitely,  and  the  system  would  be  in 
flutter.  In  the  present  situation,  we  expect  the  k-th  mode  to  have  a 
much  larger  amplitude  than  the  j-th  mode,  and  we  have  no  indication  of 
the  other  amplitudes.  Thus,  using  the  method  of  multiple  scales,  we 
obtain  a uniformly  valid  approximation  of  the  solution  of  equations 

(5.1)  in  the  form 

-6  . 

Un  * E X(T-T.’T‘>  * e ”Un1(T«*T> "T=> 

* t2S«U„!(T.,T,.TI)  ] ♦ 


• • • 


(5.2a) 
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where  the  6n  and  60  are  more  constants  to  be  determined,  and 
Tn  = en<Sot 


(5.2b) 


The  derivatives  become 

d_ 
dt 


• Do  + eSi  + e25oD2  + ... 


and 


• D20  + 2e<SoD0D1  + e2tS°(D?  + 2D0D2)  + .., 


(5.3a) 


(5.3b) 


where 


°" E ^ 


(5.3c) 

Substituting  equations  (5.2)  and  (5.3)  into  equations  (5.1)  yields 

* '*•»»'„.  + + 20.0,U„,] 

+ + <&„,  * * DiU„)+  2D,D,Uni] 

♦ 4*  [«  k V ♦ « “V*  k V 

6r  6r+6o 

+ e C2CnD0Un  ] + e L [2Cn(D0Un_+  DxUn  )] 


n'  u n,  1 n, 


* 24'CKn(D2Un<  * 0,Uni]  - [A  ^ fiJU^  * cS*Umj 


* e2{°Um)eVVS  * C j,  fnm(Umj  + c5°V  - e2S«U „.)  « 

WV 


mil 


00  00 


"1(2cosxTt)  • J,  p?,  & r"""cV  *eS,u". 


+ e26°Um  ) X (U  + e6#Un  + e2l5(,U  ) x (U„  + e6°U_ 

m2'  Po  Pi  P2  do  Pi 

+ e25oU„  ) x eNnmpq*0,  n = 1,2,...  (5.4a) 
Hi 
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where 

k = j + 1 


(5.4b) 


= “k 


(5.4c) 


and 

1 when  n = k 

■ <S-4d> 

0 otherwise 

Equation  (5.4a)  is  the  basis  for  the  analysis  of  syst  ns  having  a 
repeated  frequency,  and  we  refer  to  it  again  in  subsequent  chapters. 

In  the  present  chapter,  we  consider  the  case  in  which  there  is 
neither  a parametric  nor  an  internal  resonance.  Thus,  here  we  consider 
the  flutter  condition  for  the  linear  system  and  obtain  the  limit-cycle 
behavior  of  the  nonlinear  system^  In  this  chapter,  we  need  to  use  only  • 
a two-term  expansion,  but  in  Chapter  6 where  various  parametric  resonances 
are  studied,  we  need  a three-term  expansion. 

For  a reference,  we  use  the  j-th  mode  and  hence  let 


5j  = 0 (5.5) 

In  order  to  have  the  effect  of  damping  Included  In  a first  approx- 
imation, we  must  let 


6C  * So  (5.6a) 

And  to  account  for  the  repeated  frequency  here,  we  must  let 

5k  - Sj  3 So  or  simply  5k  * So  (5.6b) 

We  presume  that  the  j-th  and  k-th  modes  are  the  ones  primarily 
Involved  in  the  representation  of  the  flutter  system.  Thus  we  focus 
on  these  mode. 
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It  appears  that  to  bring  the  effect  of  the  aerodynamic  loading  into 
the  first  approximation,  we  must  put 

6.  + 6 - 6_  3 So 

A n m 

We  note  that,  when  n=j  and  m*k,  this  statement  leads  to 

SA  « 260  (5.7) 

and  when  n=k  and  m=j,  it  leads  to 

5.-0 

A 

If  SA  = 0,  the  aerodynamic  term  in  the  equation  governing  the  Uj  is 
elevated  to  0(e  °),  which  is  unacceptable.  On  the  other  hand,  if 

5a  = 2 S0,  such  unacceptable  result  follows;  thus,  we  let  6,  = 260. 

To  account  for  the  nonlinear  terns  in  the  first  approximation,  we 
must  put 

6N  + 6n  ' 6m  ‘ 6p  ' 6q  = 60 

When  n=j  and  m=p=q=k,  this  leads  to 

6n  = 46 o (5.8) 

For  other  combinations  of  modes,  one  obtains  different  results  for  6N, 
but  these  all  lead  to  inconsistencies. 

Here  we  do  not  consider  the  effect  of  a parametric  excitation. 

Thus,  for  the  time  being,  we  consider 


f „ = 0 for  all  n and  m 
mn 


(5.9) 
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Finally,  when  n / j or  k,  we  let 

0<  6„<S.  (5.10) 

5.1  The  Case  of  Modes  Having  Repeated  Frequencies 

A.  The  Linear  Response  - The  Flutter  Condition 

In  the  absence  of  both  an  internal  and  a parametric  resonance 
the  only  modes  remaining  in  the  analysis  are  the  j-th  and  k-th  modes. 
For  the  linear  response,  we  let 


rnmpq  = 0 for  a11  n*  m’  P- 

(5.11) 

Substituting  equations  (5.5)  thru  (5.11)  into  equation 

(5.4)  and  equating 

coefficients  of  en<5°,  we  obtain 

DoU-o  + (jD?Ui0=  0 

J J J 

(5.12a) 

DoUko  + ooj  Uk  o = 0 

(5.12b) 

DoUj  i + to  J U j x = -2DoDiUjO  ■ 2CjDolljO  + AS  j k^j^  o 

(5.12c) 

Dollki  + (A>J Uk i = - 2D0D1Uko  - Ujo  ■ 2CkDoUko 

(5.1 2d ) 

The  solution  of  equations  (5.12a)  and  (5.12b)  can  be  written  In  the  form 
Ujo  = A j ( T i ) exp(i<DjT0)  + cc  (5.13a) 


Uk°  = Ak(T1)exp(iu)jT0)  + cc  (5.13b) 

Substituting  equations  (5.13)  into  equations  (5.12c)  and  (5. 1 2d ) , and 
separating  the  result  into  real  and  imaginary  parts,  we  obtain  the 
following  conditions  for  the  elimination  of  secular  terms 


1 
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- * W * A®jkAk  ■ 0 

(5.14a) 

- 2b(V + W - aj  ■ 0 

(5.14b) 

Writing  and  Ak  as 

Aj  = Bjexp(XTi)  and  Afc  = Bkexp(XTi) 

(5.15) 

and  substituting  into  equation  (5.15)  into  equations  (5.14)  leads  to 
2icuj ( X + Cj)Bj  - = 0 (5.16a) 

8.'  + 2iajj  (X  + Ck)Bk  = 0 (5.16b) 

Setting  the  determinant  of  the  coefficients  equal  to  zero,  we  obtain 

1 21u>j<X  + Ck) 

or 

A8ilf 

X2  + (C.  + C.  )X  + (C.C.  - ) = 0 (5.17b) 

J k J k v 

Hence 

C.  + Ck  , / ABik 

»■- -ljAt  1 7(cj-ck>'  (5-17c) 

From  (5.17c)  we  see  that  the  solution  decays  if 

ABik 

<cj+ck)2>  (cj  • ck)z  +lvr  {5-18a) 

J 


= 0 , (5.17a) 


J 
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t 


or 


A > 4wjCjCk/6jk  (5.18b) 

is  the  condition  for  flutter.  Thus  in  the  absence  of  damping,  any  value 
of  A will  cause  flutter.  Here  we  should  note  that  Bjk  is  a positive 
aerodynamic  coefficient  that  is  associated  with  the  so-called  "piston- 
theory"  approximation.  Then  A is  a detuning  parameter  that  is  associated 
with  the  flutter  speed;  A is  zero  when  the  airspeed  equals  the  flutter 
speed  for  the  undamped  system.  When  A is  greater  than  zero,  the  airspeed 
is  above  the  flutter  speed.  A more  detailed  explanation  will  be  given  in 
the  numerical  example  of  Chapter  9. 


B.  The  Nonlinear  Response  - Limit  Cycle  Behavior 

To  examine  the  effect  of  the  nonlinearity  on  the  res'ponse,  we 
no  longer  use  equation  (5.11).  Again  substituting  equations  (5.5)  thru 
(5.10)  into  equation  (5.4a)  and  equating  coefficients  of  e , we  obtain 

Doll,  o + w?U . o = 0 (5.19a) 

J J J 

D0Uk„  + a3?Uko  = 0 (5.19b) 

D§Uj!  + WjUj i = - 2D0 D i U j o - 2CjDoUj0  + rjkkkuk<>  + AejkUk° 

(5.19c) 

D§Ukl  + ujUki  • - 2D0D1Uko  - 2CkD0Uk0  + Uj0  (5.19d) 

The  solution  of  equations  (5.19a)  and  (5.19b)  is  eouations  (5.13). 
Substituting  equations  (5.13)  into  equations  (5.19c)  and  (5.1 9d ) , 
Introducing  the  polar  form  (3.1),  and  separating  the  result  into  real 
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and  imaginary  parts,  we  obtain  the  following  conditions  for  the  elim- 
ination of  secular  terms: 


•3  AfJ  . . 

“j(aj  + cJaJ>  + rjkkkak  * ~T~  ak)s'"Y  - 0 

(5.20a) 

“jlak,Vi>,laisi"’'  = 0 

(5.20b) 

3 ^8  • , 

“jaj“j  * % rjkkkak  * ~r~  ak)c0SY * 0 

(5.20c) 

aiak“j  - 1 ajC0SY  * 0 

(5. 20d) 

where 


Y = cij  - ak  (5.21a) 

For  steady-state  solutions  aj  = a^  = 0.  We  see  immediately  that  a^  = ak 
= 0 are  possible  solutions.  For  nontrivial  solutions,  a^  ^ ak  t 0,  it 
follows  that  y1  = 0 then 


a'.  = (5.21b) 

The  equations  governing  the  steady-state  response  then  become 


3 ^®i  k 

“jcjaj + (s  rjkkkak + ~r~  ak)siny  = 0 

(5.22a) 

3 ^®ik 

ajaj“j  + % rjkkkak  + “2  ak^C0SY  = 0 

(5.22b) 

a .siny 

U)jCkak  + ~1 = 0 

(5.22c) 

a. cosy 

akakuj  - “'1 0 

(5.22d) 

We  now  multiply  (5.22a)  by  Ckak  and  (5.22c)  by  C^aj,  subtract,  and 
obtain 
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, A6iW  C.a? 

[(5  rjkkkak  + ~2”  )Ckak  ’ L ]sinY  = 0 

Similarly  we  multiply  (5.22b)  by  ak  and  (5.22d)  by  aj, 
obtain 


(5.23a) 
subtract,  and 


rjkkka 


2 

k 


(5.23b) 


The  coefficient  of  siny  in  equation  (5.23a)  and  the  coefficient  of  cosy 
in  equation  (5.23b)  cannot  be  zero  simultaneously.  Moreover,  siny  and 
cosy  cannot  be  zero  simultaneously.  It  follows  from  equation  (5.22c) 


that 


siny  t 0 (5.24a) 

Thus,  we  must  have 

cosy  = 0 (5.24b) 

and 

$ rjkkkak  + ASjk^Ckak  " Cjaj  = 0 (5.24c) 

If  cosy  = 0,  then  from  (5.22c)  it  follows  that 

siny  = - 1 (5.24d) 

Substituting  (5.24d)  into  (5.22c)  leads  to 

aj  * 2u)j C^a^  (5.25a) 

Then  substituting  (5.25a)  into  (5.24c)  leads  to 
MCiCk  ~ Aejk 
3rjkkk 


(5.25b) 


! 
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If  Fjkkk  is  negative  (this  Is  the  case  in  the  numerical  example 
considered  in  Chapter  9),  the  condition  for  the  existence  of  a real 
solution  is 


A > 4u)jCjCk/S.k  (5.25c) 

in  full  agreement  with  equation  (5.18b),  which  predicts  the  onset  of 
flutter.  The  linear  analysis  merely  predicts  growth,  while  the  nonlinear 
analysis  predicts  the  development  of  a finite-amplitude  lim't  cycle. 

The  existence  of  limit  cycles  was  found  numerically,  but  nrt  analytically, 
by  Dowell  (1966,1970).  Smith  and  Morino  (1976)  developed  limit  cycles 
analytically  but  not  in  the  presence  of  parametric  excitation. 


5.2  Summary 

The  analysis  of  linear  and  nonlinear  systems  having  one  repeated 
frequency  has  been  carried  out  in  the  absence  of  either  an  Internal  or 
a parametric  resonance.  The  amplitudes  of  all  the  modes  not  associated 
with  the  repeated  frequency  are  shown  to  decay.  For  the  two  modes 
associated  with  the  repeated  frequency,  the  linear  analysis  predicts  an 
infinitely  growing  solution  If  the  condition 

A > 4<ujCjCk/Sjk  (5.26) 

is  met.  The  nonlinear  analysis  requires  that  the  same  condition  be  met 
if  nontrivial  solutions  are  to  exist.  However,  the  nonlinearity  puts  a 
bound  on  the  response  and  limit-cycle  behavior  Is  found. 

Now  we  can  write  a complete  generalization  of  the  solvability 


condition  of  equation  (2.9)  that  is  valid  for  the  cases  where  either 
distinct  or  repeated  frequencies  exist.  It  has  the  form 
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2i 


i oo  ( D i A + C A ) + A l ynnAX  + I„  + Er>  + Kn  = 0 
n n n n n p"^  'np  p p n n n 

(5.27a) 


where 


and 


'JP 


pi,  AVp 

if  n + j.  m 

pi,  aVp  * Ak 

pf-1  YJPAP^P 

if  n = j 

30 

00 

pi,  aVp  * AJ 

y v.  ax  - a. 
pS]  kp  P P J 

if  n = k 

3rjkkk 

if  p = k 

_2^rjkpp  + rjpkp  + rjppk 

if  P i k 

(5.27b) 


(5.27c) 


and 


rkp  = 


3rkjjj 


L2^rkjpp  + rkpjp  + rkppj 


if  p = j 
if  P t j 


(5.27d) 


It  should  be  noted  that  for  the  case  of  distinct  frequencies  the  Kn  term 


is  set  equal  to  zero  and  equation  (5.31a)  becomes  identically  equation 
(2.9).  When  repeated  frequencies  are  considered  and  appropriate  scaling 
parameters,  6^,  are  specified  by  equations  (5.5)  through  (5.11), 
then  we  note  that  in  equation  (5.27a)  the  term 


A„  I y„.aX 

n np  p p 


P=1 


is  a higher-order  term  and  hence  does  not  enter  the  analysis  at  this 
level  of  approximation.  In  addition  equation  (5.27b)  simplifies  to 


CHAPTER  SIX 


PARAMETRIC  RESONANCE  IN  THE  ABSENCE  OF 
INTERNAL  RESONANCE 

In  this  chapter  we  investigate  several  cases  of  parametric  resonance 
in  the  absence  of  an  internal  resonance.  We  demonstrate  that  using  a 
two-term  expansion  used  in  the  method  of  multiple  scales  is  inadequate 
and  a three-term  expansion  is  required  for  the  combination  resonance. 

We  recall  from  Chapter  5 that  the  modes  not  involved  in  an  internal 
or  a parametric  resonance  decay  as  t -*■  <».  Here  the  emphasis  is  on 
steady-state  solution;  thus,  we  do  not  consider  these  modes  further. 

6.1  The  Case  Where  X ~ 2w. . 

The.  resonant  condition  associated  with  this  case  is 
X = 2ui.  + e^p  (6.1) 

J 

where  6 is  a constant  to  be  determined  and  p is  a detuning  parameter. 

Based  on  equations  (5.27a)  and  (5.27e),  we  see  that  the  solvability 
conditions  for  the  elimination  of  secular  terms  have  the  form 

- 21wj (Aj  + CjAj)  + Ij  + Ej  + 3TjkkkA2kAk  + A3j kAk  (6.2a) 

- 21cdj (Ak  + CkAk)  + Ik  + Ek  - Aj  - 0 (6.2b) 

where  I and  En  are  terms  resulting  from  internal  and  parametric  resonances, 
respectively.  In  the  last  chapter,  we  analyzed  these  equations  with  Ip  ■ 

En  * 0.  To  examine  the  effect  of  the  above  resonant  combination,  we 
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50 

must  include  the  |Cu|co  term  from  e<luat'ion  (5-4a)  with  n = j.  We  keep 
= 0 and  Ij  = Ik  = 0.  The  condition  for  inclusion  of  this  parametric 


resonant  term  is 

6F  + - 6k  = <50  (6.3a) 

Substituting  equations  (5.5)  and  (5.6b)  into  equation  (6.3a)  we  have 
5p  = 260  (6.3b) 

Hence  the  parametric  resonant  term  becomes 

Ej  = - fjk^kexp[(ie6op)T0]  (6.4a) 

Because  the  interaction  occurs  in  the  second  terms  of  the  expansion,  we 
put 

6 = 60  (6.4b) 

so  that  using  equation  (5.2b)  we  obtain 

Ej  = - fjk\exp(ipTi)  (6.4c) 


The  solvability  conditions  for  the  elimination  of  secular  terms  become 

- 12“j(Aj  * cjAj>  - fJk\«p(1PT>>  * 3rjkkkAi\  * AVk  * 0 

(6.5a) 

- f2u)j(Ak  + CkAk)  - Aj  * 0 (6.5b) 

Introducing  the  polar  form  (3.1)  and  separating  the  result  into  real  and 
Imaginary  parts,  we  obtain 

uj(aj  + Cjaj}  + aks1n8  + (l  rjkkkak  + A-^Jsiny  * 0 

(6.6a) 


l 
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“jajaj  ■ ~i~  akcos8  + (E  rjkkkak  + A-V-)C0SY  * 0 

(6.6b) 

ai 

«j(at  + Ckak)  + if-  siny  3 0 

(6.6c) 

“jVk  - r cor<  * 0 

( 6 . 6d ) 

where 

6 = pTi  - ak  - aj 

(6.7a) 

Y = aj  " ak 

(6.7b) 

For  steady-state  conditions,  aj  = ak  * 0.  Either  both  aj 

zero,  or  neither  is  zero.  For  nontrivial  solutions. 

and  ak  are 

6'  = P - aj  - - 0 

(6.8a) 

and  * 

y ■ cj  - <4  - o 

(6.8b) 

which  yield 

aj  = ak  = £ 

(6.9) 

The  steady-state  equations  become 

“jcjaj ♦ # v1n»  ♦ <1  rjkkkak  * )s1nY  ■ 0 

(6.10a) 

“jckak  * r s1n*  ■ 0 

(6.10b) 

“jaj  f - -4*  voss  + <1  rjkkkak  * A >cos*  ■ 

> 0 

(6.10c) 

wjakp  “ ajC0SY  * 0 

(6.10d) 
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It  is  possible  to  combine  equations  (6.10b)  and  (6.1  Od ) to  obtain 
aj  - wja2  (p2  + 4C2  ) (6.11a) 

Equations  (6.10b)  and  (6.1 Od ) may  be  substituted  into  (6.10a)  and 
(6.10c),  respectively,  so  that,  upon  solving  for  sin6  and  cosS  and 
squaring  the  results,  we  can  use  (6.11a)  to  obtain  a closed-form  solution: 


3T^7  Kk  * “3<4cjck  - pa> 1 -Jfs k - Spl(cj  * ck>!] 

(6.11b) 


We  note  that,  if  a nontrivial  solution  exists  as  f.^ 
hence  equation  (6.11b)  reduces  to 


-*•  0,  then 


P 


0 and 


(6.11c) 


a^  * ZuijC^a^  (6. lid) 

in  full  agreement  with  equation  (5.25),  which  gives  the  amplitudes  of 
the  limit  cycle  for  the  unforced  response. 


More  discussion  of  the  results  can  be  found  in  Chapter  9 where  a 
numerical  example  is  presented. 

6.2  The  Case  Where  X ~ ay  + aij. 

The  resonant  combination  associated  with  this  case  is 
X * uy  * Uj  + e^p  (6.12) 


' 
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To  determine  the  effect  of  the  resonant  combination,  we  must  include  the 
following  terms  in  (5.4a):  fjrUro  when  n = j and  f^U^o  when  n * r.  We 
now  have  three  equations  to  analyze  since  the  r^-mode  has  been  drawn 
into  the  response  by  the  resonant  combination  of  eouation  (6.12).  If  the 
scaling  conditions  specified  in  equations  (5.5)  thru  (5.10)  are  met 
along  with  equation  (6.3b),  we  see  that  the  resonant  term  of  the  form 
fjnUr°  in  ecluat10n  (5.4a),  for  n = j,  will  be  of  higher  order  and  will 
not  enter  the  two-term  expansion.  Satisfying  equation  (5.10)  for 

3 6o  allows  for  inclusion  of  the  f-  II  0 term;  however  the  term 
r jr  r 

frkUr°  in  e<luat10n  (5.4a),  for  n = r,  will  be  of  higher  order  and  will 
not  enter  the  two-term  expansion.  In  this  case  the  solution  predicts 
that  Ap  decays  as  t > ».  Thus  it  appears  that  the  original  two-term 
expansion  is  an  inadequate  approximation,  and  recourse  Is  made  to  a 
three-term  expansion. 

We  re-examine  equation  (5.4a)  and  proceed  with  an  analysis  which 
parallels  that  of  Chapter  5. 

x 

We  arbitrarily  choose  the  scaling  coefficients,  e °,  so  that  the 
effects  of  damping,  the  repeated  frequency,  the  aerodynamic  loading,  the 
nonlinearity,  and  the  parametric  resonance  all  interact  at  the  same 
order  of  approximation.  We  have  already  determined  that  this  will  not 

x 

occur  at  0(e  °),  so  we  expect  that  a uniformly  valid  solution  will  exist 
at  0(e*5°).  Thus  the  6 In  equation  (6.12)  is  set  at 

6 * 26a  (6.13a) 


54 


For  a reference,  we  again  use  the  j-th  mode  and  hence  let 

6-  = 0 (6.13b) 

J 

In  order  to  include  the  effect  of  damping,  now  we  must  let 


5c  ■ 2S0 


(6.13c) 


Accounting  for  the  effect  of  the  repeated  frequency,  we  must  let 


sk  - 5j  ■ 24» 


Using  (6.13a)  we  obtain 


5.,  = 250 


(6.13d) 


To  bring  the  aerodynamic  loading  into  the  approximation,  we  set 

4A  *4n-5m=24«  . 

where  we  note  that  for  n = j and  m = k 

6a  = 46 o (6.13e) 

and  for  n = k and  m = j 

4a’° 

Following  the  reasoning  of  Chapter  5,  we  conclude  that  equation  (6.13e) 
must  be  satisfied  to  produce  acceptable  results. 

To  account  for  the  nonlinearity,  we  set 

6W  + 6 - 6 - 6 - 6 * 260 
n n m p q 


When  n 3 j and  m * p ■ q * k,  this  leads  to 
SN 


6„  * 86 o 


(6.1 3f ) 
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Other  combinations  of  n,  m,  p,  and  q lead  to  different  values  for  6N, 
but  all  these  cases  lead  to  inconsistencies  in  the  expansion. 

Finally,  to  include  the  effect  of  the  parametric  resonance,  we  set 

6F  + 6n  ' 5m  = 260 
Setting  t 

6p  = 36  0 (6. 13g) 

allows  for  acceptable  solutions.  If  n = j and  m = r,  then  equation 
(6.13g)  yields 

Sr  = 60  (6.1 3h ) 

and  one  of  the  parametric  resonance  terms  f . U 0 can  be  accomodated. 

jr  r 

When  n * r and  m = k the  other  resonant  terms  frkUko  is  also  Included. 
Substituting  equations  (6.13)  into  equation  (5.4a)  and  equating  co- 


efficients of  equal  powers  of  en<So,  we  obtain 

D§U.„  + (d?U -o  = 0 (6.14a) 

J J J 

D§Uko  + tojUko  = 0 (6.14b) 

D§Uro  + wJ;Uro  - 0 (6.14c) 

Dollji  + coj  Uj  1 + 2D0DiUj0  + 2fjkUk0cosXT0  = 0 (6.15a) 

Dollki  + oj?Uki  + 2D0DiUkoa  0 (6.15b) 

0$Uri  + o)^Uri  + 2D0D1Uro  3 0 (6.15c) 
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D§U,z  + u)?U.2  + 2D0D2U,o  + DfU.0  + 20oDiU  • 1 + 2C-D0U  • , 

JJJ  J J J JJ 


* 2fj,.Ur.cosXT,  * 2fJkUklcosXT,  - rJkkkUj. 


' ^jkV  ■ 0 

D§Uk2  + U)jUk2  + 2D0D2Uko  + D?Uko  + 2D0D1U|<i  + 2CkD0Uk0 
+ Ujo  - 0 

0§Ur2  + a£Ur2  + 2D  o D2Ur  o+  D?Uro  + 2000^  + 2CrD0Uro 

+ 2frkUkocosXT0  = 0 

The  solutions  of  equations  (6.14)  are 
Uj°  = Aj(T1,T2)exp(iu)jT0)  + cc 

Uk°  = Ak(Ti,T2)exp(ia)jT0)  + cc 

Uro  ■ Ar(T1,T2)exp(iu)rT0)  + cc 

Substituting  equations  (6.17)  into  equations  (6. 15), we  find 
terms  are  eliminated  provided 


DiAj  = DtAk  * - 0 

and  that  the  An  are  independent  of  Ti. 
It  follows  that 


V * Tjk 


Akexp[i(x  + (Dj  )T0]  7fkexp[i(x  - ajjT0] 

* + x(x  - '2a ,T 

J 


X(X  + 2uij) 


(6.16a) 

(6.16b) 

(6.16c) 

(6.17a) 
(6.17b) 
(6.17c) 
that  secular 

(6.18) 

+ cc 


(6.19a) 
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and 

Ukl  « Url  = 0 (6.19b) 

Substituting  (6.17)  through  (6.19)  into  (6.16)  yields  the  solvability 
conditions 

- 2i(VAj + W • fjAexp(ipT2)  + 3rjkkkAk\  + A6jkAk = 0 

(6.20a) 

- 2iojj (A^  + CkAk)  - Aj  = 0 (6.20b) 

- 2iur(A^  + CrAr)  - frk\exp(iPT2)  = 0 (6.20c) 

where  primes  denote  differentiation  with  respect  to  T2.  Letting 

Aj  = \ an(T2)expfiotn^T2^  n * J*  k»  r (6.21 ) 

ft 

substituting  equation  (6.21)  into  equations  (6.20)  and  separating  the 
result  into  real  and  imaginary  parts,  we  obtain 


'ir  ^ A6*. 

“j ( a j + Cjaj)  + + arsinSl  + l*  rjkkkak  + + ak>si"Y  " 0 
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These  equations  do  not  readily  admit  a closed-form  solution.  They 
are  discussed  further  in  Chapter  9,  where  a numerical  example  is  pre- 
sented. We  note  that  only  five  equations  are  independent,  equation 
(6.25f)  being  a combination  of  equations  (6.25e)  and  (6.25d).  A sixth 
equation  can  be  obtained  from  equations  (6.23). 

S2  * Si  + Y (6.25g) 

6.3  The  Case  Where  X ~ mr  - a-. 

In  this  case  the  ordering  scheme  turns  out  to  be  the  same  as  that 
of  the  previous  section.  Hence,  we  write  the  resonant  combination  as 


X = (Jr  - coj  + e2^°p  (6.26) 

The  resonant  terms  become 

EjexpOuijTo)  = fj^expCKuj  - e2%)T0]  (6.27a) 

and 

Erexp(1u)rT„)  - frkAkexp[i Ur  + e2%)T0]  (6.27b) 

with 

01  8 - pT2  + ar  - (6.28a) 

Y ■ aj  - ak  (6.28b) 

02  8 pT2  - op  + ak  (6.28c) 


The  governing  equations  have  the  same  form  as  those  of  the  last  section. 
Hence  the  steady-state  equations  become 
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^ir  3 

ujCjaj  + “T  arsinSl  + rjkkkak  + “T"  ak]siny  * 0 (6-29a) 


“jcjai  * r slBY  ■ 0 


u C ] + a, 

r r r 2 k 


rk 


a.sInSz 


fr ah 


o " cos62  + rrfT~  cosy 


r r 


5^rk 


(6.29b) 

(6.29c) 

(6.29d) 


° ■ 2^  cosSl  * 3J5J  { % arcosfli  - [|  rJkkk«J 


A6, 


* -r  ak]cos^ 


(6.29e) 


a.cosy  1 I fir  3 A0ik  1 

* ^aj  1 -r  v°s6'  - c s rjkkkak  + -r  akJ“^< 


(6. 29f ) 

6.4  Summary 

Several  cases  in  which  parametric  resonances  are  involved  have  been 


investigated.  The  resonance  case  X ^ 2w,  appears  at  first  order,  while 

J 


the  resonance  cases  X ^ ay  ± Uj  appear  at  second  order.  Nothing  con- 
clusive is  self-evident  from  a cursury  examination  of  the  steady-state 
equations.  Analysis  and  discussion  await  the  coverage  of  numerical  examples 
in  Chapter  9. 


CHAPTER  SEVEN 

THE  EFFECT  OF  INTERNAL  RESONANCE 


Again  we  choose  to  investigate  an  internal  resonance  of  the  form 
us  T 3oj..  In  this  chapter  we  examine  parametric  resonances  in  the  presence 

' J 

of  this  internal  resonance.  We  also  discuss  the  case  where  the  internal 
resonance  is  present  but  there  Is  no  parametric  resonance,  the  only 
excitation  coming  from  the  aerodynamic  loading. 

7.1  The  Case  Where  X Z 2us.  and  u»„  ~ 3us. . 

J r J 


The  resonances  associated  with  this  case  are 


X = 2uj.  + c p 

J 


and 


u)„  = 3us.  + e u 
r J 


(7.1a) 


(7.1b) 


In  this  case,  the  general  form  of  the  solvability  conditions  is  as 
follows: 

’ 2iwj(Aj  + CjV  + + Ej  + 3rjkkkVk  + A8jkAk  = 0 (7-2a) 


- 21u)j (A£  + CkAk)  + Ik  + Ek  - Aj  = 0 (7.2b) 

- 21«r(A;  + CrAr)  + Ir  + Er  - 0 (7.2c) 

To  obtain  the  best  approximation,  we  need  to  develop  an  ordering 
scheme  which  includes  the  effect  of  the  repeated  frequency,  the  parametric 
resonance,  the  internal  resonance,  and  as  many  other  terms  as  it  Is 
consistent  to  retain. 
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Proceeding  as  in  Section  6.1,  we  find  that  equations  (5.5)  thru 
(5.10)  and  (6.3b)  still  hold  and  that 

«r  s 0 (7.3a) 

Thus,  the  nonzero  resonant  terms  in  the  solvability  conditions  become 


Ej = - fjk\exp<l(,T1> 

(7.3b) 

Er  * ■ frkAkexp^p  “ a^Tl^ 

(7.3c) 

and 

'r  ' rrkkkAkexp(-,oT'> 

(7.3d) 

Substituting  equations  (7.3b)  - (7.3d)  into  equations 

(7.2),  introducing 

the  polar  form  (3.1),  and  separating  the  result  into 

real  and  imaginary 

•parts,  we  obtain 

7 AS-. 

wj(aj  + cjaj}  + “V"  sinBl  + (s  rjkkkak+  ~T 

ak)siny  - 0 

(7.4a) 

fikak  3 ASik 

- -4^  c°s«>  * rjkkkak  * ~r  ak>«^  * 0 *7,4bi 

ai 

°jj(ak  + Ckak)  + -f-  siny  = 0 

(7.4c) 

ai 

“jak“k  - r * 0 

(7.4d) 

“r(ar  * crar>  ' ff  rrkkkak  s,n“  * lfrkaks,n8!’  0 

(7.4e) 

VA  + S rrkkkakcos®!“  T^-kV056*  * 0 

where 

(7. 4f ) 

Bi  a pTi  - ak  - aj 

(7.4g) 

Y * aj  - ak 

(7.4h) 
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V = - aTi  + 3ak  - ar 


62  = (p  - a)T!  + ak  - ar  = Si  + Y + U 


(7.41) 

(7.4J) 


For  the  steady-state  solutions  aj  = a£  a aj.  = 0,  and  it  follows  that 

lution* 

(7.5a) 


a-  = a,  = a =0  are  possible  solutions.  For  non-trivial  solutions 
j k r 


Bi  * 02  3 y'  = Y1  » o 


so  that 


and 


0 * ak  * “j 


a'.  = a! 
J k 


0 = 3ak  ar 


from  which 

“j  = »k  ■ f 


and 


ar  = - 


(7.5b) 

(7.5c) 

(7.5d) 

(7.5e) 

(7.5f) 


The  steady-state  equations  become 
f..,a 


s1"6'  + (f  rjkkkak  * -T— >s1n*  = 0 (7-6a) 

a. 

wjCkak  + 2 siny  = 0 (7.6b) 

i frkaks1n82 

»rCrar  - J WS*1"*  * ' 0 (7'6c) 

_ ^ilcak  1 ^®ikak 

wjai  ‘ 2 cos^  + rjkkkak + -j—^y  - 0 


(7.6d) 
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w.a.p  - a. cosy  * 0 (7.6e) 

J K J 

i , frkakCOS62 

Vr(f  -'JV 0 <7*6f> 

Equations  (7.6a,  b,  d,  e)  are  identical  to  equations  (6.10)  and  hence 
can  be  reduced  to 

ak  * 3fT^  C-A8jk  + - p*)  ± /fjk  - fejp'tCj  * Ck)!  ] 

(7.6g) 

aj  = cojak  /p2  + 4Ck2  (7.6h) 

Moreover,  8i  and  y can  be  obtained  from  these  equations.  To  find  af,  62 
and  u,  we  solve  equations  (7.6c,f)  and  ( 7 . 4 j ) numerically.  The  response 
is  the  same  as  that  in  Chapter  6 for  the  two  modes  associated  with  the 
repeated  frequency.  Due  to  the  presence  of  an  internal  resonance,  a 
third  mode  is  drawn  into  the  response,  but  it  does  not  influence  the 
j-th  and  k-th  modes. 

7.2  The  Case  Where  A Z m + u.  and  ^ 3io,. 

L J L 1 

The  resonances  associated  with  this  case  are 


X 


+ “j  + 


e5p 


(7.7a) 


and 


, . 6 

+ £ o 

r J 

The  solvability  conditions  of  equations  (7.2)  remain  valid, 
with  the  formulation  of  the  three-term  expansion  of  Section 


(7.7b) 

Starting 
6.2,  we  note 
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that  the  parametric  resonances  have  been  accounted  for  and  that  we  must 
attempt  to  include  the  internal  resonance  terms 

‘j  * (rjrkk  * rjkrk  * rjkkr>Vk  <7'8a> 

and 

!r  = rrkkkAk  exp("ie<SaTo)  (7.8b) 

2 x 

Unfortunately,  both  of  these  terms  are  order  e 0 and  do  not  enter  at 
this  level  of  approximation. 

We  thus  conclude  that  the  effects  of  the  parametric  resonance  and 
the  repeated  frequency  dominate  the  response  and  that  the  effect  of  the 
internal  resonance  is  of  higher  order.  To  this  order  of  approximation 
the  response  in  the  presence  of  a parametric  resonance  is  the  same 
regardless  of  whether  an  internal  resonance  is  present  or  not. 

7.3  The  Case  of  No  Parametric  Resonance  with  ui  ~ 3aij. 
r J 

Only  an  internal  resonance  is  associated  with  this  case  and  It  is 
associated  with  the  combination 

ay  = 3u>j  + e^a  (7.9) 

The  solvability  conditions  of  equations  (7.2)  remain  valid.  There  Is  no 
parametric  resonance  so  the  En  3 0.  Using  the  scaling  conditions  of 
equations  (5.5)  thru  (5.10),  (6.3b)  and  (7.3a),  we  find  that  the  only 
Internal  resonance  term  that  appears  to  first  order  is 

'r  • rrkkkAk  «PHc5oT.) 


(7.10) 


66 


r 


Substituting  equation  (7.10)  into  equation  (7.2),  introducing  the  polar 
form  (3.1),  and  separating  the  result  into  real  and  imaginary  parts,  we 
obtain  equations  identical  to  equations  (7.4)  when 

(7.11) 


fjk  :frks0 


The  resulting  steady-state  equations  lead  to  the  following: 


4VCjCk  - A6jk 


3rjkkk 


(7.12a) 


aj  = 2ujCkak 


(7.12b) 


and 


1 r 

ft  * 


8 rkkJTk 


.A’fc/  * 7) 


(7.12c) 


Hence  we  note  that  equations  (7.12a)  and  (7.12b)  are  the  limit-cycle 
equations  developed  in  Chapter  5 when  no  resonance  was  present.  The 
modes  associated  with  the  repeated  frequency  are  independent  of  the 
internal  resonance,  which  only  serves  to  draw  the  r-th  mode  into  the 
response.  Moreover,  the  steady-state  amplitudes  a^  and  ak  predicted  by 
equations  (7.6g)  and  (7.6h)  approach  those  given  by  equations  (7.12a) 
and  (7.12b)  as  f^k  vanishes.  We  note  that.  If  a steady-state  solution 
exists,  |fjk|  > | p | 2oj j ( C j + Ck)  and  hence  p must  vanish  as  fjk  vanishes. 

Finally,  we  note  that  equations  (7.6c)  and  (7.6f)  yield  equation  (7.12c) 


as  fjk  and  fpk  vanish. 
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7.4  Summary 

We  have  examined  two  cases  of  parametric  resonance  in  the  presence 

of  the  internal  resonance  u>  ~ Zu-. 

r J 

a.  When  X is  near  twice  the  repeated  frequency  the  response  for 
the  repeated  modes  is  the  same  as  for  the  case  of  no  internal  resonance. 
The  third  mode  is  merely  drawn  into  the  response. 

b.  When  X is  away  from  twice  the  repeated  frequency  the  response 
tends  to  be  the  same  regardless  of  the  presence  of  an  internal  resonance. 
The  effect  of  the  parametric  resonance  and  the  repeated  frequency  dominate 
the  response. 

c.  In  the  case  of  no  parametric  resonance  or  a parametric  resonance 
away  from  twice  the  repeated  frequency,  the  form  of  the  response  for  the 
modes  associated  with  the  repeated  frequency  is  the  same  as  that  discussed 
in  Chapter  5 where  no  resonance  was  present.  Hence  the  only  mode  affected 
is  the  mode  drawn  into  the  response  by  the  internal  resonance. 


CHAPTER  EIGHT 


LATERAL  VIBRATIONS  OF  COLUMNS 


The  lateral  vibrations  of  columns  are  used  as  an  example  where  the 
analysis  contained  in  Chapters,  2,  3,  and  4 can  be  applied.  The  natural 
frequencies  of  oscillations  are  distinct.  The  governing  equation  will 
be  nondimens ional ized  and  converted  to  the  form  discussed  in  Chapters  2, 
3,  and  4. 


8.1  General 


Relatively  large-amplitude  lateral  deflections  cause  significant 
stretching  of  the  neutral  axis  when  the  longitudinal  displacements  of 
the  ends  are  restrained.  For  small,  but  finite  deflections,  one  needs 
to  use  nonlinear  geometrical  relationships  to  account  for  this  stretching. 


The  equations  governing  the  deflection  can  be  written  in  the  form 

£ 

9w 

at 

(8.1) 


E.f^pAt^  ■[*.)♦!*/  <$>*«£  -C 


The  cubic  term  accounts  for  the  stretching.  The  terms  which  account  for 
shear  deformation  and  rotary  inertia  are  formally  the  same  order  as  the 
term  which  accounts  for  stretching.  However,  these  terms  are  linear, 
and  consequently,  they  are  responsible  for  only  small  perturbations  of 
the  solution  of  the  linear  problem.  Here  the  focus  is  on  nonlinear 
effects . 


It  is  convenient  to  Introduce  nondimensional  variables  (denoted  by 
primes);  thus,  we  let 
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1 


I 


r2 

x = Lx 1 , w = — w , 

t = F-Vf  *' 

(8.2a) 

p(t)  = EA(ff  p'(t'), 

2Ar3  r~r 

c ~rz~  VpE 

C' , e=  -pr  EA 

(8.2b) 

Substituting  equation  (8.2)  into  equation  (8.1)  and  dropping  the  primes, 
we  obtain 


(8.3) 


Here  the  characteristic  length  is  one-half  the  actual  length  of  the 
beam.  We  seek  approximate  solutions  of  equation  (8.3)  which  are  valid 
for  small,  but  finite,  values  of  e. 

We  express  the  deflection  as  an  expansion  in  terms  of  the  linear 
free-osci llation  modes  as 

w(x,t)  « l Um(t;e)4>m(x)  (8.4) 

m=l 

In  equation  (8.4),  the  <|>m  are  the  orthonormal  solutions  of  the  following: 

4>miv  - ud$_  =0  (8.5a) 

m mm 

where 

<t>m  * 0 and  ^ = 0 (8.5b) 

at  clamped  ends  and 

$ = 0 and  ^ = 0 (8.5c) 

at  hinged  ends.  The  are  the  natural  frequencies  and  the  4>m  are  the 
linear  free-osci llation  modes. 

Substituting  equation  (8.4)  into  equation  (8.3),  multiplying  by  <f>n, 


and  integrating  the  result  from  x s 0 to  x = 2,  we  obtain 


f 
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00  00  QO 


Jn  + <*£U„  = e[-2C0SXt  7 fU  - 2C„  Cl„  + y y T r UUU] 

n n n m=i  nm  m n n ^ ^ A nmpq  m p qJ 


for  n = 1 , 2,  ... 


(8.6a) 


where 


P(t)  = F cosXt 

2 2 

2f  = - F I <h*"dx  = F f dx 

nm  J n m J nrm 


and 


nmpq 


4.  C 

/ W*  / 


*p»qdx 


(8.6b) 


(8.6c) 


Equations  (8.6a)  and  (2.1)  are  the  same. 
8. 2 Hinged-Hinged  Beam 
In  this  case. 


4>n  = sin  x 


(8.7a) 


and 


^ = l n2Tr2 


(8.7b) 

It  follows  that  resonant  combinations  of  natural  frequencies  abound; 
however,  for  this  beam 
1 


T s * t u u 6 6 

nmpq  4 n p nm  pq 


(8.8a) 


Thus  all  the  I are  zero.  Essentially,  there  is  no  internal  resonance 
for  this  beam. 


From  equation  (8.6b)  it  follows  that 

f _ * i w F6 
nm  Z n nm 


(8.8b) 
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Therefore,  resonances  of  the  combination  type,  X near  o)n  ± u)m,  are  not 
possible  for  this  beam. 

Using  equations  (8.7)  and  considering  X to  be  near  2oon , we  obtain 
the  following  from  equations  (3.7): 

an  + ^nan  + 1 Fansiny  = 0 (8.9a) 

anan  " 3ajH/32  ' 1 FanC0Su  = 0 (8.9b) 

where 


U = pTi  - 2otn 

For  the  non-trivial  steady-state  solution,  we  find  that 

4C  16p  - 3(jj„a^ 

siny  = — p—  , cosy  - gp 

and 


It  follows  from  (8.10b)  that 


F>  4Cn 

if  a real  solution  for  an  can  exist.  If  F _>  4Cn,  then 
( 1 ) when 

0 < - \ /f"2  - 16C*  , 

*vo  real  solution  exists; 

2 ) 

- I A1  - 16C*  <o£7  ^ p2  - 16Cn 


(8.9c) 


(8.10a) 


(8.10b) 


(8.11) 


(8.12a) 


(8.12b) 
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one  real  solution  exists; 

(3)  when 

p > | / F2  - 16C*  (8.12c) 

two  real  solutions  exist. 

New  we  consider  the  stability  of  the  various  possibilities  . Using 
the  analysis  in  Chapter  3,  equations  (3.8)  through  (3.13),  we  find  that 
for  this  case  equation  (3.11)  becomes 

n = - C ± /Cl  ± Jp-  a*  / F2  - 4C*  (8.13) 

n ^7  n 32  ns  n 

When  the  negative  sign  inside  the  radical  is  used,  the  real  parts  of 
both  values  of  are  negative  and  hence  the  coresponding  motion  is 
stable,  the  negative  sign  corresponds  to  the  larger  ap  given  by  equation* 
(8.10b).  Hence,  when  two  real  solutions  exist,  only  the  one  with  the 
large  amplitude  is  stable,  and  when  only  one  real  solution  exists,  it  is 
stable. 

The  stability  of  the  trivial  solutions  can  be  examined  in  a similar 
fashion.  It  turns  out  that,  when  p is  in  the  region  where  there  is  no 
real  solution  of  equation  (8.10b),  the  trivial  solution  is  stable;  when 
p is  in  the  region  where  there  is  one  real  solution  of  equation  (8.10b), 
the  trivial  solution  is  unstable;  and  when  p is  In  the  region  where 
there  are  two  real  solutions  of  equation  (8.10b),  the  trivial  solution 
is  stable. 

The  results  are  illustrated  in  Figures  1-3  when  \ is  near  ui.  For 
combinations  of  F and  p in  Region  I,  no  nontrivial  steady-state  solution 
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P 

Figure  1.  Parameter  space,  F is  the  amplitude  and  p(ep  * X - 2u>i)  is 


the  frequency  of  the  excitation. 
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exists,  and  the  trivial  solution  is  stable.  In  Region  II,  the  trivial 
solution  is  unstable,  while  the  only  nontrivial  solution  is  stable. 

In  Region  III,  two  nontrivial  solutions  are  possible.  The  trivial 
solution  and  the  larger-ampl itude  nontrivial  solution  are  stable,  while 
the  smaller-ampl itude  nontrivial  solution  is  unstable. 

In  Region  I,  it  appears  that  regardless  of  how  large  the  initial 
amplitude  is,  the  phasing  never  becomes  such  that  the  rate  at  which  work 
is  being  done  is  as  large  as  the  rate  at  which  energy  is  being  dissipated. 
Indeed,  if  the  damping  term  were  deleted  from  equation  (8.9a),  the 
predicted  response  would  still  decay.  Thus,  we  conclude  that,  for  the 
linear  as  well  as  the  nonlinear  system,  the  phasing  is  such  that  the 
force  actually  does  negative  work  and  contributes  to  the  decay. 

In  Region  II,  for  a linear  system  the  phasing,  which  does  not 
change  with  amplitude,  is  such  that  work  is  being  done  at  a faster  rate 
than  energy  is  being  dissipated  and  thus  the  response  to  any  initial 
disturbance  grows  without  bound.  For  the  nonlinear  system,  the  phasing 
for  large  amplitudes  differs  from  the  phasing  in  the  linear  system  as  a 
consequence  of  the  nonlinear  term  in  equation  (8.9a).  The  effect  is  to 
limit  the  rate  at  which  work  is  being  done  to  the  rate  at  which  energy 
is  being  dissipated  and  thereby  to  produce  a bounded  harmonic  response. 

If  the  initial  disturbance  is  very  large,  the  response  will  decay  until 
the  steady-state  solution  is  reached,  while  if  the  initial  disturbance 
is  very  small  the  response  will  grow  (the  system  being  governed  by  the 
linear  equations  when  the  amplitude  is  small)  until  the  nonlinear  term 
in  equation  (8.9b)  becomes  large  enough  to  cause  the  phase  to  shift 
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significantly.  Thus,  in  Region  II  all  initial  disturbances  produce  the 
same  steady-state  response. 

In  Region  III,  the  response  of  the  linear  system  to  an  initial 
disturbance  always  decays.  The  mechanism  causing  the  decay  is  the  same 
as  in  Region  I.  The  results  for  the  nonlinear  system  show  that  only  the 
larger  of  the  two  possible  steady-state  responses  is  stable.  Thus,  it 
appears  that,  if  the  initial  disturbance  is  small,  the  nonlinear  term  in 
equation  (8.9b)  does  not  have  a strong  influence  on  the  resulting 
motion  and  the  system  behaves  essentially  as  a linear  system;  the  motion 
decays.  On  the  other  hand,  when  the  initial  disturbance  is  large  enough 
to  produce  the  nontrivial  steady-state  response,  the  nonlinear  term  has 
a strong  influence,  and  changes  in  the  phase  occur.  Thus,  in  Region  III 
there  is  the  possibility  of  producing  motions  which  have  characteristics 
that  are  similar  to  those  of  the  motions  in  Region  I as  well  as  those  in 
Region  II.  The  type  of  motion  is  determined  by  the  amplitude  and  phase 
of  the  initial  disturbance.  This  is  a rare  example  in  which  a nontrivial, 
steady-state  response  of  the  nonlinear  system  exists  in  a region  where 
the  response  of  the  linear  system  decays. 

Let  us  suppose  that  the  frequency  of  the  excitation  is  varied  while 
the  amplitude  is  held  fixed.  This  process  is  represented  by  the  line 
through  points  Aj,  A2  and  A3.  Between  points  A!  and  A2,  the  trivial 
solution  is  stable;  moreover,  it  is  the  only  steady-state  solution 
possible.  Between  points  A2  and  A3,  the  trivial  solution  is  unstable, 
and  the  only  realizable  solution  is  the  one  given  by  equation  (8.10b). 
Beyond  point  A3,  the  trivial  solution  is  again  stable  and  two  solutions 
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are  realizable  --  the  trivial  solution  and  the  larger  of  the  two  solutions 
given  by  equation  (8.10b). 

Finally,  let  us  suppose  the  amplitude  of  the  excitation  Is  varied 
while  the  frequency  is  held  fixed.  This  process  is  represented  by  the 
line  through  points  B3,  B2  and  B3.  Between  points  B3  and  B2,  only  the 
trivial  solution  is  possible  because  the  requirement  of  equation  (8.11) 
is  not  satisfied.  Between  points  B2  and  B3,  two  solutions  are  realizable-- 
the  trivial  solution,  which  is  stable,  and  the  larger  of  the  two  solutions 
given  by  equation  (8.10b).  Beyond  point  B3,  the  trivial  solution  is 
unstable  and  the  only  realizable  solution  is  the  one  given  by  equation 
(8.10b). 

For  the  first  process,  the  amplitude  of  the  solution  is  plotted  as 
a function  of  p in  Figure  2.  If  the  frequency  of  the  excitation  decreases 
from  a large  value,  we  note  that  upon  reaching  point  A3,  where  the 
trivial  solution  becomes  unstable,  there  is  a jump  in  the  value  of  ai  up 
from  point  A3  to  point  A3.  This  process  is  indicated  by  the  arrows. 

For  the  second  process,  the  amplitude  of  the  solution  Is  plotted  as 
a function  of  F in  Figure  3.  We  note  the  presence  of  a jump  phenomenon. 

If  F is  increased  slowly  from  Point  0 toward  point  B3,  a3  remains  zero. 
However,  the  trivial  solution  becomes  unstable  at  B3  and  a slight 
increase  in  F at  this  point  causes  ai  to  jump  up  to  point  Bj.  Then 
further  increases  cause  a3  to  follow  curve  CB3D  toward  point  D.  When  F 
is  decreased  a3  follows  curve  CB3D  away  from  point  D pasts  point  Bj, 
without  jumping  down  to  point  B3,  until  it  reaches  point  C.  At  this 
point  a slight  decrease  in  F causes  ai  to  jump  back  to  zero  at  point  B2. 


REGION  II 
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The  arrows  indicate  this  path.  Figure  3 corresponds  to  p greater  than 
zero.  As  p decreases,  point  C approaches  point  B3,  and  the  multi-valued 
region  and  hence  the  jump  vanish. 

The  results  for  A near  2 o)n,  where  n = 2,  3 are  similar  to 

the  above  results. 

Next  we  consider  an  example  for  which  the  internal  resonance  plays 
an  important  role. 


8.3  Hinged-Clamped  Beam 

* 

For  a beam  which  is  hinged  at  the  left  end  and  clamped  at  the  right 

sin2/“ojI 


*n  = 


n fsin  /uT 
n I n 


x - 


n 


sinh  vuT 


sinh2vruT 


n 


n x 


(8.16a) 


where 


sinh2^  oT 


K = 


n i / 

[sinh22/  uj  - sin22/~oT  ] /2 


(8.16b) 


and  the  o>„  are  the  roots  of 
n 

tanh2>4on  - tan2/  u)n  = 0 (8.16c) 

The  first  five  roots  of  equation  (8.16c)  are 

cui  = 3.855,  W2  = 12.491,  0)3  = 26.062,  u)i»  s 44.567,  and 

0)5  = 68.008 

Referring  to  equation  (2.8a),  we  find  that  for  the  first  five  modes  the 
combinations  for  which  ea  is  below  eight  percent  of  o>n  are 

3o>i  s 0)2  - eoi  (8.18a) 


I 


0)  i ■*"  0>2  CO  3 = 0)4  ■*"  CO  2 


2u>2  * Wi*  = 0)5  ■*■  CO  3 


oj2  + 2 0)3  = 0)5  + eow 


0 It,  + 0)3  ■ 0)1  = 0)5  + £Os 


2o)4  “ 0)3  = 0)5  + £06 


(8.18b) 
(8.18c) 
(8. 18d) 
(8. 18e) 
(8.1 8f ) 


Referring  to  equations  (2.9),  we  find  that  the  conditions  for  the 
elimination  of  secular  terms  are 

0)1  (a  1 + C 1 a i ) - 3Qia ia2sinui  + Q2a2a3a4Siny2  ~ Qsa3a4asSiny5 


- Imag[Eiexp(-iai )]  = 0 


(8.19) 


oiiaial  + 4-  ai  7 y + 3Qiaia2cosui  + Q2a2a3a4Cosy2 
8 im  m 


+ Q5a3a4a5COSu5  + Re[Eiexp(-iai )]  = 0 


(8.19b) 


oi2(a2  + C2a2)  + Qia^sinyi  + Q2a!a3a4Sinu2  + 2Q3a2a4a5siny3 


+ Q4a3a5siny4  - Imag[E2exp(-ia2)]  = 0 


(8.20a) 


oj2a2a2  + i a2  7 y a*  + Qiafcosyi  + Q2aia3a4COsy2 
B m=1  2m  m 

+ 2Q3a2a4a5siny3  + Q4a3a5siny4  + RE[E2exp(-1a2)]  * 0 

(8.20b) 

o>3(a3  + C3a3)  + Q2aia2a4siny2  + 2Q4a2a3a5s1ny4  + Q5a1a4a5siny5 


- Q6aia5siny6  - Imag[E3exp(-ia3)]  ■ 0 


(8.21a) 
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w3a3oi3  + g a3  l y ma*  + Qza^a^cosya  + 2Q4a2a3a5cosy4 
m=l  3 

+ Q5aia4a5cosus  + Q6aJasCOSy6  + Re[E3exp(-ioi3)]  = 0 

(8.21b) 

ujulai  + C4a4)  - Q 2a i a2a 3 si n]j2  + Q3a2asSinu3  + QsajajasSinys 

+ 2Q6a3a4a5sinu6  - Imag[E4exp(-ici4)]  = 0 (8.22a) 

s 

oj4a4ai  + i-  a4  V y a*  + Q2a:a->a3cosy2  + Q3a2ascosy3 

O ___  i u 111  111 

m=  I 

+ Q5aia3a5cosy5  + 2Q6a3a4a5cosy6  + Re[E4exp(-ia4)]  3 0 

(8.22b) 

us(as  + C5a5)  - Q3a2a4Siny3  - Q4a2a3Siny4  - Q5aia3a4siny5 

- Q6a3a4Siny6  - Imag[E5exp(-1a5)]  = 0 (8.23a) 

5 

wsascis  + ^ a5  l y a*  + Q3ala4Cosy3  + Q6a2a3Siny4 
m=l 

+ Q5aia3a4COSys  + Qeasa^cosye  + Re[E5exp(-ia5)]  = 0 

(8.23b) 

where 

Qi  3 r i i 12/8 

Q2  3 ( r 1423  + r1  2 3 4 + 3 24 )/4 

Q3  3 (r2254  + 2r2425)/8 

Q4  3 (r2  5 3 3 + 2r  2 3 3 5 )/8 

Os  3 (r  14  3 5 + T1345  + H 5 34  )/4 


Q 6 3 ( r 3 5 4 4 + 2r*  3 4 4 s ) /8 
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Ui  = ct2  - 3ai  + aiTi 

U2  = Qt  1 “ Cl2  ” 01 3 - 0l4  “ C2Ti 

y3  = Za2  + a4  - as  + a3Ti 

U 4 = 02  + 2a3  - as  + O^i 

U5  = ct3  + Cti*  - Ol  “ 05  “ O5T1 

y 6 = 2a4  - ot3  - as  + OsTi 
and  the  En  depend  on  the  excitation. 

In  this  section,  the  following  three  cases  are  considered: 

(1 ) X near  2on 

(2)  X near  2oj2 

(3)  X near  an  + an 

It  follows  that  the  Ep,  for  n » 3,  4,  and  5,  are  zero  for  all  three 
cases  being  considered;  thus,  from  equations  (8.19)  - (8.23),  it  follows 
that  a possible  solution  is 

a3  = a4  = a5  = 0 

and  a3  and  a2  are  given  by  „ 

o>i (ai  + - 3Q1a?a2siny1  - Imag[Eiexp(-iai)]  = 0 (8.24a) 

onaiol  + ^ ai(Ynaf  + mal)  + 3Qia?a2cosy!  + ReCE^xpHa! )]  = 0 

(8.24b) 

and 

0)2 (a 2 + C2a2)  + QiaiSlnyj  - Imag[E2exp(-ia2)]  = 0 (8.25a) 

o)2a2a2  + ^ a2(Y2iai  +Y22a2)  + Qiajcosyi  + Re[E2exp(-1o2)]  * 0 

(8.25b) 
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which  are  either  equations  (4.30)  or  (4.45)  depending  on  the  form  of  the 
parametric  resonance.  Tso  and  Asmis  considered  the  system  above,  but 
took  Qj  to  be  zero. 

Case  (1),  A near  2m 

To  express  the  nearness  of  X to  2u>i,  we  introduce  another  detuning 
factor  p defined  by 

X = 2u)i  + ep  = 0)2  - wi  + e(p  - ai)  (8.26) 

where  equation  (8.18a)  was  used.  Then  equations  (8.24)  and  (8.25)  can 
be  rewritten  exactly  as  (4.30)  with  m=l  and  k=2. 

an  (a  i + C^i)  + Jr  fuajsing  - ^ fi2a2sin(S  - ui) 

- 3Q!aia2sinui  3 0 (8.27a) 

oo i a iot  1 + g-  ai(yi j a i + yizal)  - g-  fnaicoss  - g-  fi2a2C0s(S  - ui) 


+ 3Q!aia2COSui  = 0 (8.27b) 

o)2(a2  + C2a2)  + g-  f 2 iaisin(6  - ui)  + Qia^sinui  = 0 (8.27c) 

o)2a2a2  + g-  a2(y2iai  + y22ai)  - g-  f2iaiCOs(s  - Ui)  + Qiaicosui  3 0 


where 


(8.27d) 

6 = pTi  - 2a  i 

(8.28a) 

Ui  = aiTi  + a2  - 3a i 

(8.28b) 

For  the  steady-state  response, 

al  3 a2  3 0 and  p{  3 8'  3 0 

(8.29) 
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It  follows  that  the  trivial  solution  is  possible: 

a i = a 2 = 0 (8. 30 ) 

If  neither  ai  nor  a2  is  zero,  one  can  use  equations  (8.28)  and  (8.29)  to 
rewrite  equations  (8.27)  as 

cj i C i a i f i iaiSin8  ■ fi232sin(8  - p2)  - 3Qiaia2sinpi  = 0 

(8.31a) 

coiaxp  + -g-  (yi ia?  + Yi232)ai  - ^ fuaicosg  - ^ fi2a2cos(;3  - ui) 

+ 3Qiaia2cosui  = 0 (8.31b) 

w2C2a2  + ^ f2la!Sin(S  - p2)  + Qialsinpi  = 0 (8.31c) 

3 1 1 

u>2az(j  p - oi)  + g (Y22af  +.  Y22ai)a2  - 2 f2ja2cos(6  - y3 ) 

+ Qiaicospi  = 0 (8. 31 d ) 

A Newton-Raphson  technique  is  used  to  find  the  solutions  of  equations 
(8.31). 

In  Figure  4,  some  typical  results  show  the  amplitudes  of  the  first 
and  second  modes  (ai  and  a2)  as  functions  of  the  amplitude  of  the  ex- 
citation F,  (refer  to  equation  8.6b).  We  note  the  multiplicity  of  possible 
jumps;  jumps  are  indicated  by  the  arrows.  We  can  trace  the  histories  of 
a2  and  a2  as  F increases  very  slowly  from  zero  to  a large  value  and  then 
slowly  decreases  back  to  zero.  Initially,  both  ai  and  a2  are  zero,  and 
they  remain  zero  until  point  A is  reached.  At  this  point  the  trivial 
solution  becomes  unstable.  Then  the  amplitudes  jump  either  to  curve 
BCDE  or  FGH.  If  the  jump  is  to  G,  then  as  F continues  to  increase,  the 
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for  constant  frequency  of  the  excitation,  p(ep 


response  follows  curve  FGH.  If  the  jump  is  to  D,  then  as  F continues  to 
increase  the  response  follows  BCDE  until  another  jump  occurs  from  E to 
H.  An  additional  increase  in  F causes  the  response  to  follow  FGH  again. 
Decreasing  F at  this  time  causes  the  amplitudes  to  remain  on  FGH  until  F 
is  reached.  Then  a jump  to  either  C or  the  trivial  solution  occurs.  If 
the  jump  is  to  C,  then  an  increase  or  decrease  in  F causes  the  response 
to  follow  BCDE  until  the  next  jump  is  encountered.  The  jump  at  E has 
already  been  discussed  and  the  one  at  B returns  the  amplitudes  to  the 
trivial  solution.  In  Figure  5 the  amplitudes  are  shown  as  a function  of 
the  detuning  p.  Possible  jumps  are  again  indicated  by  the  arrows.  The 
response  for  a!  is  initially  similar  to  that  of  Figure  2 beginning  with 
the  trivial  solution  when  p is  very  negative.  As  p increases  the  trivial 
solution  becomes  unstable  and  the  amplitudes  follow  curve  ABC  as  in 
Figure  2.  As  p increases  beyond  C a transfer  of  energy  occurs  between 
the  modes,  and  a2  grows  at  the  expense  of  ai  until  E is  reached.  The 
dominance  of  the  second  mode  is  quite  clear.  An  additional  increase  In 
p produces  jumps.  The  amplitudes  jump  either  to  G or  they  jump  to  the 
trivial  solution.  A further  increase  in  p will  cause  the  response  to 
continue  along  FG  or  remain  trivial,  depending  on  which  way  the  jump 
occured. 

At  this  point,  if  the  trivial  response  was  attained,  the  solution 
would  remain  trivial  for  a decrease  in  p until  H is  reached.  Here 
another  jump  occurs  from  H to  B since  the  trivial  solution  Is  unstable, 
below  H.  Once  at  B,  an  Increase  or  decrease  In  p would  cause  the  response 
to  follow  ABCDE. 


p(ep  = X - 2u)i ),  for  constant  amplitude  of  the  excitation. 
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If  the  jump  from  E had  been  to  G and  then  p decreased,  the  response 
would  follow  GF  to  F where  jumps  again  occur.  Both  amplitudes  could 
become  trivial  or  jump  to  D.  An  increase  or  decrease  in  p at  this  point 
would  cause  the  response  to  follow  ABODE  again  or  remain  trivial  until  H 
is  reached  as  discussed  above. 

Case  (2),  X near  2oj2 


In  this  case,  we  introduce  the  detuning  parameter  p according  to 


X = 2oj2  + Ep 

(8.32) 

For  the  steady-state,  equations  (8.34)  and  (8.35) 

reduce  to  equations 

(4.49)  with  m = 1 and  k = 2 

3i(ojiCi  - 3Q ja i a2si ny i ) 3 0 

(8.33a) 

ai[o>i(ai  + pV3  + -g  (yii^i  +Yi2a2)  + 3Qia 

ia2cosui]  3 0 

(8.33b) 

u)2C2a2  + f22a2sine  + Qiaisinui  3 0 

(8.34a) 

<j2pa2  + (y2iai  + y22a2)a2  - y f22a2cos8 

+ Qiaicosui  3 0 

where 

(8.34b) 

8 = pTi  - 2a2 

(8.35a) 

and  as  above 

Ui  * OiT ! + a2  - 3at 

(8.35b) 

From  equations  (8.33)  and  (8.34),  It  appears  that  there  are  two 
possibilities  for  a nontrivial  solution:  either  a!  is  zero  and  a2  is 
nonzero  or  both  ai  and  a2  are  nonzero.  Thus,  the  situation  Is  the  same 
as  that  In  Case  (1).  When  at  is  zero,  the  results  resemble  those 
illustrated  in  Figures  1-3;  consequently,  they  are  not  given  here. 
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The  results  when  both  ai  and  a2  differ  from  zero  are  illustrated  in 
Figures  6 and  7.  The  top  and  bottom  branches  go  together,  and  the  middle 
two  branches  go  together.  The  important  thing  to  note  is  that  over  a 
fairly  wide  range  a!  is  much  larger  than  a2.  This  means  that  the  first 
mode  dominates  the  deflection,  though  the  second  mode  is  the  only  one 
directly  excited.  Consequently,  when  both  a!  and  a2  differ  from  zero, 
there  can  be  a significant  transfer  of  energy  from  the  second  mode  to 
the  first. 

Case  (3),  X near  uii  + m> 

In  this  case  we  introduce  a detuning  parameter  p defined  as 

X = 4wi  + ep  = w2  + tui  + e(p  - a i)  (8.36) 

Also,  we  let  # 

Pi  = aiTi  + a2  - 3ai  (8.37a) 

U2  = (p  - ai)Ti  - on-  a2  (8.37b) 

For  steady-state  solutions  a^  = v'n  - 0,  and  equations  (8.24)  and  (8.25) 
yield  those  of  (4.22)  with  m = 1 and  k = 2;  that  is 

Co i C x a i + ^ 2 2 s1ny2  - 3Qia2a2sinyi  = 0 (8.38a) 

c*»iax  f + g-(Yua?  + YizaDaj  - cosy2  + SQ^azCOSy!  = 0 

(8.38b) 

oo2C2a2  + siny2  + Qiafslnyi  * 0 (8.39a) 

(^a^!2-  - oj)  + g-  (Yi2a?  + Y22a|)a2  - cosu2  + Q1aJcosm  * 0 

(8.39b) 
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Two  possibilities  exist:  either  ai  and  a2  are  zero,  or  neither  Is 
zero.  Figures  8 and  9 show  the  variation  of  ai  and  a2  with  F and  p.  In 
both  figures,  the  top  curves  for  ai  go  with  the  bottom  curves  for  a2, 
and  the  two  middle  curves  go  together.  Again  we  note  that  over  a wide 
range,  aj  is  much  larger  than  a2,  and  hence  the  first  mode  dominates  the 
response. 

8.4  Observations  - Internal  and  Parametric  Resonances. 

a.  The  Effect  of  the  Detuning  Parameter  p on  the  Response, 
an  vs  F,  X near  2ui. 

Some  interesting  observations  may  be  noted  concerning  the  behavior 
represented  in  Figure  4,  which  depicts  the  relationship  between  the 

response  amplitudes,  a , and  the  amplitude  of  the  forcing  function,  F. 

n • 

The  detuning  parameter  p measures  the  closeness  of  the  frequency  of  the 
excitation  to  2u>i.  In  Figure  10,  when  p is  small,  the  response  is 
essentially  indistinguishable  from  that  of  Figure  3,  i.e.  the  case  when 
no  internal  resonance  is  present.  The  amplitude  of  the  second-mode  is 
almost  zero.  As  the  frequency  of  the  excitation  increases,  (i.e.,  p 
increases),  the  effect  of  the  internal  resonance  begins  to  appear.  The 
response  may  resemble  that  in  Figure  10,  or  the  second  mode  may  dominate 
the  first  mode  as  shown  in  Figures  11  and  12.  As  the  detuning  of  the 
parametric  resonance  approaches  that  of  the  internal  detuning  (p  = o), 
we  see  In  Figure  11,  strong  modal  Interaction  and  the  possible  dominance 
of  a2  over  for  a large  range  of  the  amplitude  of  the  excitation,  F. 

In  Figures  4 and  12,  for  p greater  than  a,  the  second  can  dominate  the 
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Figure  8.  The  modal  amplitudes  as  functions  of  the  amplitude 
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first  mode,  over  a very  limited  range  of  the  amplitude  of  the  ex- 
citation. However,  strong  modal  coupling  is  present,  and  the  second 
mode  amplitude  becomes  relatively  significant  at  the  expense  of  the 
first;  hence  we  note  energy  being  transferred  from  the  excited  mode  to 
the  one  involved  in  the  internal  resonance.  Figures  4,  10,  11,  and  12 
graphically  demonstrate  the  transition  from  a single-mode  response  to 
a multimode  response.  The  larger  the  detuning  p is,  the  more  the  second 
mode  enters  the  response. 

b.  The  Effect  of  Damping  and  the  Amplitudes  of  the  Excitation 
on  the  Response,  an  vs  p. 

Using  Figure  5 as  a point  of  departure,  we  see  that  the  effects  of 
damping  on  the  system  can  be  made  apparent.  As  a general  rule,  the 
smaller  the  damping  coefficients  Cn  a're,  the  more  interaction  one  can 
expect.  Increasing  the  damping  coefficients  tends  to  reduce  the  modal 
interaction  to  the  point  where  the  second-mode  response  is  below  that  of 
the  first  mode,  and  there  is  no  longer  a crossover  point  between  a x 
and  a2.  In  addition,  jump  phenomena  tend  to  disappear  as  damping 
increases.  Figures  13  and  14  demonstrate  both  of  these  effects.  As  the 
damping  increases,  the  response  approaches  tnat  of  Figure  2 as  a limit, 
or  the  response  aproaches  the  single-mode  response  without  internal 
resonance.  Increasing  the  amplitude  of  the  excitation  tends  to  spread 
the  stable  and  unstable  branches  and  straighten  the  curves,  first  the 
unstable  branch,  then  the  stable  one.  We  note  this  behavior  in  com- 
paring Figure  13  with  Figure  14  and  Figure  5 with  Figure  15.  As  F 
Increases  beyond  the  zone  where  jumps  can  occur,  (the  central  portion 
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Figure  13.  The  modal  amplitudes  as  functions  of  the  frequency 
of  the  excitation,  p(ep  = A - 2a)i),  for  constant 

amplitude  of  the  excitation,  F.  stable, 

unstable. 
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of  Figure  4),  then  a further  increase  in  F tends  to  force  a single-mode 
response,  approaching  that  of  Figure  2. 

8.5  Summary 

The  method  of  solution  described  in  Chapters  2,  3,  and  4 is  used  to 
solve  two  numerical  examples  — one  without  internal  resonances  (a 
hinged-hinged  column)  and  one  with  internal  resonances  (a  hinged-clamped 
column). 

For  hinged-hinged  columns,  there  are  in  effect  no  internal  resonances 
and  no  combination  resonances.  Thus,  the  parametric  excitation  can,  in 
the  first  approximation,  excite  only  one  mode.  Jump  phenomena  can  be 
produced  by  varying  either  the  frequency  or  the  amplitude  of  the  excita- 
tion. In  this  case  the  frequency  of  the  response  is  always  one-half, 
that  of  the  excitation. 

For  hinged-clamped  columns,  there  are  both  internal  resonances  and 
combination  resonances.  When  X is  near  2o>i  or  on  + W2,  both  the  first 
and  second  modes  are  strongly  excited  (i.e.,  appear  in  the  first  approx- 
imation), but  generally  the  first  mode  dominates  the  second.  Qualitatively, 
the  response  in  both  cases  closely  resembles  a single-mode  (the  first) 
response.  However,  in  the  first  case  (X  near  2a)!)  the  frequency  of 
the  response  is  one-half  that  of  the  excitation,  while  In  the  second 
case  it  is  one-fourth  the  frequency  of  the  excitation. 

When  X is  near  2o)2,  two  nontrivial  solutions  are  possible  — 
either  the  second  mode  only  is  strongly  excited  or  the  second  and  first 
modes  are  strongly  excited.  In  the  second  case,  the  first  mode  domi- 
nates, and  the  response  qualitatively  resembles  the  response  when  X is 
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near  2uh  and  when  X is  near  oji  + u2.  Multiple  jumps  are  possible.  When 
only  the  second  mode  is  excited,  the  frequency  of  the  response  is  one- 
half  that  of  the  excitation;  on  the  other  hand,  when  both  modes  are 
strongly  excited,  the  frequency  of  the  response  is  one-sixth  that  of  the 
excitation. 

Thus,  for  the  cases  considered,  it  is  possible  to  produce  responses 
having  frequencies  that  are  one-half,  one-fourth,  and  one-sixth  that  of 
the  excitation. 

One  excitation  can  produce  rather  dissimilar  responses.  On  the 
other  hand,  several  different  excitations  can  produce  responses  which 
appear  to  be  quite  similar. 

When  X is  near  2un»  detuning  the  excitation  allows  the  internal 
resonance  to  induce  strong  modal  coupling  and  the  second  mode  can  be 
made  to  dominate  the  response.  Increasing  the  damping  and/or  the 
amplitude  of  the  excitation  inhibits  modal  coupling.  If  the  increase 
is  large  enough,  it  can  result  in  essentially  a single-mode  response. 


CHAPTER  NINE 


PANEL  FLUTTER 


The  analysis  contained  in  Chapters  5,  6,  and  7.  is  applied  to 
the  motion  of  a plate  in  an  airstream  undergoing  cylindrical  bending  (no 
spanwise  bending).  As  in  the  last  chapter,  the  governing  equations 
are  nondimensionalized,  but  here  they  are  converted  to  the  form  discussed 
in  Chapters  5,  6,  and  7.  The  process  is  slightly  more  involved  than  for 
the  beam. 


9.1  General 


Following  Dowell 
as  follows 

n 3“w  . h a2w 

D + pm  h W 


(1966,  1975),  we  write  the  governing  equation 


(||)2dC] 


3w  P»U»  r 3w  /M2  - 2 > 1_  3w  , 

3t  ‘ (M2  - 1)1/2  3x  M2  - 1 U»  9t  J 


(9.1) 


where  D is  the  flexural  rigidity,  pm  Is  the  material  density,  h is  the 
panel  thickness,  a is  the  panel  length,  p^  is  the  density  of  the  air  in 
the  frees tream,  is  the  frees tream  velocity,  and  M is  the  Mach  number, 
and  Nx  refers  to  the  inplane  parametric  loading  and  is  harmonic  in  what 
follows.  The  edges  of  the  plate  are  restrained  so  that  stretching 
occurs,  and  the  term  which  accounts  for  it  has  the  same  cubic  form  as 
that  of  a beam.  Modal  damping  is  again  assumed.  The  final  term  accounts 
for  the  aerodynamic  pressure  due  to  the  plate  motion  and  it  is  in  the 
form  of  the  so-called  "piston- theory"  approximation. 
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It  is  convenient  to  introduce  nondimensional  variables  (denoted  by 
primes);  thus  we  let 


1 


x = p ax' , 

x 4D 


w = 


t . a2  pmh 


6p  Nva2 
£ R - 


N = 2e\.  =_Jl_c 

a 4 /Thb 


(9.2) 


The  parameter  which  is  a measure  of  the  dynamic  pressure  is 
p U2a3 


X = 


8(M2-1)1/2D 

The  coefficient  of  the  last  aerodynamic  term  is 
(M2-2)  2*p~a 


(m2-i)9/“  pi 


T 


(9.3) 


(9.4) 


and  can  be  taken  as  small  because  the  ratio  of  the  density  of  air  to  the 
density  of  the  material  is  small. 

Substituting  equations  (9.2)  and  (9.3)  into  equation  (9.1)  and 
dropping  the  primes,  we  obtain 


0 


9 Cr  3w 
- 2e  C 9t 


(9.5) 


where  60»  <Sp»  <5^,  etc.  conform  with  the  ordering  of  Chapters  5-7. 

We  express  the  deflection  as  an  expansion  in  terms  of  the  linear  free- 
osci llation  modes;  that  is 


w(x,t)  = Um(t;e)4>m(x) 


m=l 


(9.6) 
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The  d>  are  the  orthonormal  solutions  of 
Tm 

♦m™  - ^m  = 0 

where 

6 - 0 and  = 0 (9.7a) 

m m 

at  clamped  edges  and 

0=0  and  0"  = 0 (9.7b) 

m m 

at  hinged  edges.  The  ojm  are  the  natural  frequencies  and  the  0m  are  the 
linear  free-oscil lation  modes. 

Substituting  equation  (9.6)  into  equation  (9.5),  multiplying  by  0 , 
and  integrating  the  result  from  x = 0 to  x = 2,  we  obtain 


•• 

— 00 

6 

p “ 

Un  + 

co*U  + X l 

n n 

m- 1 

Ct„mUm  = [-  £ 

nm  m L 

J 2cosXtf  Um 
nm  m 
m- 1 

6r  . 

s 1 l 

m=1  p=l 

00 

■ 2e  \Un  + 
n n 

qa'j  rnmpqUfnUpUq-* 

(9.8a) 

where 

2 

a„m  = 
nm 

j[  v;dx 

(9.8b) 

Rx  = 

FcosXt 

(9.8c) 

2 

£ 

2 

r 

2fnm= 

- FJ  ^n^mdx 

- f J 

‘ FGnm 

( 9 , 8d ) 

0 

0 

l 

2 • 

r 

nmpq 

s J Vi;dx 

J 0p0q<ix  = - G 

nmGpq 

(9.8e) 

0 0 


At  this  stage  with  the  beam,  we  had  achieved  the  desired  form  of  the 
governing  equations,  but  this  is  not  so  for  the  flutter  problem. 
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Equation  (9.8a)  may  be  written  in  matrix  form  as 

" '■'C 

AU  + BU  = e FU  + e CU  + e ‘TU 


(9.9) 


provided  we  use  a finite  number  of  terms  in  the  summation.  The  matrices 
U and  U are  column  matricies  of  length  M.  The  remainder  are  square 
matrices  of  dimension  M x M.  Then  we  note  that 


[A]  = [I] 

The  identity  matrix 

[B]  = (VJ  + X[ot  ] 

where  the  frequency  matrix  is  diagonal, 


[F]  = - 2cosAt[fnni] 


and 


^ ’ • 2 CcrJ 


[T]  - - HU]1  [Gnm][U][Gpq] 


(9.10a) 

(9.10b) 

(9.10c) 
(9.1 Od ) 

( 9 . 1 Oe ) 


The  system  of  equations  which  has  been  analyzed  in  Chapters  5 through  7 
is  of  the  form 

6 5 6 5 

IU-+  JU  = e FF’U  + e CC’U  + e AAB'U  + e NT'U  (9.11a) 

where  the  first  coefficient  matrix  is  the  identity  matrix  and 


[J] 


<*>1 

0 

0 0 . 

. . 0 

1 

u>? 

0 . . 

. . 0 

0 

0 

CO  3 . 

• • • 

• • • 

• • • 

0 

• 

• • • 

• • • 

• • 

(9.11b) 
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is  of  the  Jordan  canonical  form,  which,  in  this  case  Is  a diagonal 
matrix  with  one  off-diagonal  term  that  is  non-zero;  F',  C'  and  T1 
contain  the  same  type  of  terms  as  their  counter  parts  in  equation 
(9.9). 

The  matrix  B'  contains  aerodynamic  coefficients  which  appear  when 
the  dynamic  pressure  causes  the  stability  parameter  X to  be  near  the 
critical  value  XQ  which  causes  two  natural  frequencies  to  merge  and 
flutter  to  ensue.  A refinement  of  the  matrix  B'  is  required  before  an 
attempt  can  be  made  to  transform  equation  (9.9)  into  equation  (9.11a), 
the  Jordan  form.  Equation  (9.10b)  may  be  written 


[B]  = |V]  + X r«]  - ( X - X)[a  ] 
nJ  c nm  c nm 


= [B  ] + (X  - X )[a  ] 

c c nmJ 


(9.12a) 


(9.12b) 


where  X - X is  a measure  of  the  nearness  of  X to  the  critical  value 
c 

Xc,  and  can  be  expressed  in  terms  of  a detuning  parameter  as  follows 


= X - X. 


(9.12c) 


Then  it  is  actually  B which  must  be  converted  to  Jordan  form  when  the 

c 

value  of  the  stability  parameter  reaches  the  critical  value,  Xc,  and  the 
onset  of  flutter  occurs.  The  two  lowest  eigenfrequencies  of  Bc  become 
equal  at  this  time  The  larger  eigenfrequencies  remain  very  near  the  wn 
for  n > 3. 


It  is  now  possible  to  introduce  a similarity  transformation  which 
transforms  the  matrix  Bc  to  a Jordan  form.  (A  complete  discussion  of 
the  choice  of  this  matrix  is  contained  in  Appendix  A).  We  thus  choose  a 
matrix  Z such  that 


W 
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U = ZV  (9.13a) 

Then  the  similarity  transformation  that  leaves  the  eigenvalues  unchanged 


becomes 

r n 

- - 

o 

o 

o 

o 

CJ 

3 

Vi 

Z-1B  ZV  = JV  = 
c 

£ 

r> 

o 

o 

o 

v2 

0 0 u)|  ...  0 

• •••••• 

• 

1 

. 

• o 

1 

. VM  . 

Z-1AZV  = Z'lIZV 

= IV 

(9. 

and 


Z_1FZV  = F'V 


Z_1TZV  = T'V 


(9.13b) 


(9.13d) 

(9.13e) 


This  last  transformation  merits  some  detailed  consideration.  From 
equation  (9.13a)  we  can  write  the  transpose  as 

UT  = VTZT  (9.13f ) 


Then  using  equations  (9.10e)  and  (9.1 3e ) we  write 

Z_1TZV  = Z_1VTZTG  ZVGZV 
nm  pq 

Letting 


(9.1 3g ) 


( 9 . 1 3h ) 


we  note  that 

VTGiV  (9.131) 

Is  a quadratic  form  and  it  is  thus  a scalar.  Hence  we  may  write  equation 

(9.13g)  as 

VTGiVZ-1G  ZV 

pq 


(9.1 3 j ) 


no 


Letting 

G2  = Z_1G  Z 

pq 

we  write 

T'V  = - VTGiVG2V 


(9.13k) 


(9.131) 


Equation  (9.131)  has  the  same  cubic  form  as  before,  or 

1 I I V V V 

m p q nmpq  m p q 


where 


r1  = Gi  g2 

nmpq  xnm  zpq 


(9.14a) 


(9.14b) 


Finally,  the  second  half  of  the  aero-matrix  of  equation  (9.12b)  can  be 
transformed  into 


B'  = - Z"1 [a]Z  = 6 


Hence 


..  5p  6r  6. 

IV  + JV  = e r F 1 V + e ' V + e NT'V  + e aABV 


(9.15) 


(9.16) 


and  the  panel  flutter  problem  has  been  reduced  to  the  form  of  equations 
(5.1).  We  shall  now  discuss  the  response  for  the  various  cases  of 
parametric  and  internal  resonances  examined  in  Chapters  5-7. 

9.2  Simply  Supported  Plate. 

For  cylindrical  bending  of  a simply  supported  plate  we  may  write 


and 


b = sin  / co_  x 
rn  n 


“n  = l n2u2 


(9.17a) 


(9.17b) 


Ill 


where  the  <|>n  and  oin  are  the  free-oscl nation  modes  and  frequencies.  It 
follows  from  equation  (9.8)  that 

— — if  m + n is  odd 

n2  - m2 

0 if  m + n is  even 

f = u)  F<5 

nm  2 n nm 

nm  2 n nm 

r = - r r 

nmpq  ’ " nm  pq 

The  transformations  of  equations  (9.13)  are  applied  to  equations  (9.18) 


to  provide 

612  = 99.2  (9.19a) 

f12  = 404.9  2 (9.19b) 

f i 3 - - 3917  \ (9.19c) 

f32  = -1017  j (9.1 9d ) 

f33  = 22.74  \ (9.19e) 

r1222=  - 5729  (9.1 9f ) 

and 

r3 2 22s  “ 1.438  (9.1 9g ) 

The  frequencies  become 

coi  = cj2  = 8.1076  (9.20a) 

o)3  3 22.066  (9.20b) 


(9.18a) 

(9.18b) 
(9.18c) 
(9.1 8d ) 
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While  the  critical  value  of  the  dynamic-pressure  parameter  is 

Ic  = 42.917  (9.20c) 

where  the  wn  and  Xc  satisfy 

*iv  + \c<p'  - w2*  =0  (9.21) 

The  various  cases  of  Chapters  5-7  are  discussed,  beginning  with  the 
response  in  the  presence  of  a parametric  resonance.  The  cases  considered 
are 

(1 ) X near  2m 

(2)  X near  m ± m 
9.3  Parametric  Resonance 

a.  Case  (1)  X ~ 2m 


The  governing  equations  for  this  case  were  developed  in  Section  6.1, 
and  for  convenience  we  rewrite  equations  (6.11)  as 


a!  = — (-A8i2  + oji(4CiC2-  P2)  t /f?2  - 4u>ip2 (Ci  + C2 ) 2 ) 3 

jkkk 


(9.22a) 


di  s 0)1^2  yp2  +4C|  (9.22b) 

A general  observation  that  can  be  made  immediately  concerns  the  relative 
sizes  of  a!  and  a2.  If 

/p2  + 4Cf  < 1 (9.23a) 

Hence,  for  a2  > ai  the  detuning  must  be  very  nearly  zero  and  the  damping 
extremely  small,  otherwise  a!  > a2.  We  note  that  even  though  at  > a2 
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the  response  is  not  necessarily  dominated  by  the  first  mode.  The  modal 
amplitude  a:  multiplies  a one  (1)  whereas  a2  multiplies  e 0 and  hence 
its  effect  can  be  much  greater.  Figure  16  shows  the  variation  of  the 
modal  amplitude  with  the  amplitude  of  excitation  and  the  response  has  j 

essentially  the  same  form  as  that  in  Figure  3.  The  response  contains 
the  same  regions  where  trivial  and  nontrivial  solutions  exist  and  jump 
phenomena  are  possible.  We  note  that  as  the  amplitude  of  the  parametric 
resonance  term  increases,  0 <_  F _<  Fi,  the  only  solution  is  the  trivial 
one  and  it  is  stable.  Then  we  enter  the  region  where  two  nontrivial 
solutions  exist,  Fi<  F £ F2,  one  stable  and  one  unstable.  Finally  when 
F ^ F2  the  trivial  solution  is  unstable  and  only  one  nontrivial  solution 
exists,  which  is  stable,  and  a jump  occurs.  If  we  scan  in  the  opposite 
direction  we  start  with  a nontrivial  stable  solution  and  remain  on  that 
upper  branch  until  we  reach  the  region  where  the  only  solution  is  the 
stable  trivial  solution  and  a jump  occurs.  One  can  refer  to  Figure  3 
where  the  arrows  indicate  the  jump  phenomena  and  the  regions  where 
trivial  and  nontrivial  solutions  exist  are  shown.  The  form  of  the 
response  for  the  second  mode  is  the  same  as  the  first  except  the  ampli- 
tude is  so  small  that  it  appears  to  be  zero  by  comparison. 

The  analog  of  Figure  2,  modal  amplitude  vs  detuning  of  the  excita- 
tion frequency,  is  Figure  17.  This  figure  shows  the  modal  amplitude 
ai  vs  the  detuning  of  the  dynamic-pressure  parameter,  A.  However,  the 
form  of  the  response  is  essentially  the  same  and  exhibits  the  same 
regions  of  stability  and  jump  phenomena.  As  we  scan  from  large  negative 
A to  large  positive  A,  or  alternatively,  as  we  increase  the  dynamic 
pressure  from  below  the  flutter  speed,  we  encounter  first  the  region 

es- 1 


400  800  1200  1600 

F 

Figure  16.  First  modal  amplitude  as  a function  of  the  amplitude  of 


the  excitation,  F,  for  constant  frequency  of  the  excitation 
p(ep  * X - 2wi),  and  constant  aerodynamic  detuning,  A. 

stable,  unstable. 
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A 

Figure  17.  First  modal  amplitude  as  a function  of  aerodynamic 
detuning.  A,  for  constant  frequency  of  excitation 
p(ep  = X - 2uu)  and  constant  amplitude  of  excitation, 
F.  stable, unstable. 


116 


where  the  trivial  solution  is  the  only  solution  and  it  is  stable,  A <_Ai. 
Between  Ai  <_  A <_  A2  the  trivial  solution  is  unstable  and  the  only 
nontrivial  solution  is  stable.  Beyond  A > A2  there  is  again  a stable 
trivial  solution  and  two  nontrivial  solutions,  only  the  larger  one  being 
stable.  Again  the  form  of  the  second  mode  is  essentially  the  same 
except  the  magnitude  appears  to  be  zero  by  comparison. 

In  Figure  18  we  plot  the  modal  amplitude  vs  the  detuning  of  the 
parametric  excitation  p.  This  new  graph  has  a character  of  its  own  and 
has  no  analog  with  the  case  of  distinct  frequencies.  However,  it  still 
contains  regions  where  multiple  solutions  exist  and  points  out  where 
jumps  can  occur.  We  can  scan  the  frequency  axis  and  observe  the  fol- 
lowing. For  p _<  Pi  only  the  trivial  solution  exists  and  it  is  stable. 
Between  px  <_  p <_  p2  two  nontrivial  solutions  exist,  the  larger  being 
stable,  the  smaller  unstable.  In  addition,  the  trivial  solution  is 
stable,  and  since  we  are  scanning  from  the  left,  the  response  remains 
trivial.  Between  p2  <_  p <_  p3  only  one  nontrivial  solution  exists  and  it 
is  stable.  Concurrently,  the  trivial  solution  becomes  unstable  and  an 
upward  jump  occurs.  Arrows  indicate  the  jumps  and  paths  followed  by  aj. 
As  the  frequency  increases  between  p3  _<  p _<  p4  the  response  remains 
nontrivial,  since  that  solution  is  stable.  A nontrivial,  unstable 
solution  also  exists  as  well  as  a stable  trivial  solution.  For  p > p4 
we  see  that  the  only  solution  possible  is  the  trivial  one,  which  is 
stable,  and  a downward  jump  occurs  where  the  nontrivial  solution  ceases 
to  exist.  Scanning  from  the  right  we  note  a mirror-image  response.  The 
solution  is  first  trivial  until  p * p3.  It  jumps  to  the  nontrivial 
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Figure  18.  First  modal  amplitude  as  a function  of  frequency 
of  excitation,  p(ep  * X - 20)0 , for  constant 
amplitude  of  excitation,  F,  and  constant 
aerodynamic  detuning,  A.  stable,  unstable. 
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solution  and  remains  nontrivial  until  p = pi  at  which  time  it  jumps  back 
to  the  trivial  solution  for  further  decreases  in  p. 

When  the  damping  is  small  and  the  frequency  of  excitation  is  almost 
exactly  tuned,  or  p is  very  near  zero,  the  condition  of  equation  (9.23a) 
can  be  met.  The  second  mode  is  greater  than  the  first  mode  and  we  note 
from  Figures  19  and  20  that  the  characteristic  shape  of  the  curves  is 
the  same.  Hence,  the  amount  of  detuning  and  damping  will  determine 
which  mode  amplitude  is  larger. 

It  should  be  noted  that  for  the  plate  considered,  Ti222  is  negative 
in  equations  (9.22a)  and  the  existence  of  any  solution  other  than  the 
trivial  one  depends  upon  the  conditions 

A812  ■ wi  (4CiC2  ~ p2)  - J f 2 2 - 4toiP2  (Ci  + C2  )2  > 0 

(9.24a) 

and 

f f 2 > 4toip2  (Cx  + C2  )2  (9.24b) 

If  (9.24b)  is  not  satisfied  then  there  are  no  nontrivial  solutions.  The 
same  holds  for  (9.24a).  If  both  conditions  are  met  then  there  are  two 
nontrivial  solutions.  The  stable  one  is  associated  with  the  positive 
value  of  the  radical.  When  only  one  solution  exists,  it  is  also  associated 
with  the  positive  value  of  the  radical.  In  addition,  it  is  dependent  on 
the  magnitude  of  the  amplitude,  F,  of  the  parameteric  excitation,  the 
relative  sizes  of  the  detuning  and  the  damping,  and  the  aerodynamic 
detuning  parameter,  A.  We  note  that  in  this  case,  812  is  a positive 
constant  and  that  a positive  value  of  A refers  to  a dynamic  pressure  in 
excess  of  the  critical  flutter  speed.  The  presence  of  large  damping 
tends  to  necessitate  large  forcing  amplitudes,  large  detunings  of  the 
parametric  terms,  and  large  values  of  dynamic  pressures  to  produce 


Figure  20.  The  modal  amplitudes  as  functions  of  the  aerodynamic 
detuning,  a,  for  constant  frequency  of  excitation, 
p(ep  a A - 2aii ) and  constant  amplitude  of  excitation, 
F.  stable, unstable. 
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nontrivial  solutions.  However,  this  same  condition  make  a2  -very  small 
compared  with  oj.  On  the  other  hand,  small  damping  permits  strong 
interaction  of  the  modes  and  even  allows  a2  to  be  much  greater  than  ai 
when  detuning  is  small.  The  amplitude  of  the  parametric  excitation  and 
the  aerodynamic  detuning  merely  influence  the  amplitude  of  the  response 
but  not  the  relative  sizes  of  ai  and  a2. 

A few  special  cases  may  also  be  examined.  If  no  parametric 
resonance  is  present,  (9.24a)  has  the  form 

A3 1 2 ~ 4C i C2cuf  > 0 (9.25) 

and  there  is  either  the  trivial  solution  or  one  nontrivial  solution, 
depending  on  whether  (9.25)  is  satisfied.  Figure  21  illustrates  the 
response  for  small  damping.  When  large  damping  is  present  ai  > a2.  As 
C2  -*-0,  ai  -*•  0 and  the  only  mode  excited  is  the  second,  in  spite  of 
the  parametric  force  exciting  the  first  mode.  The  response  is  as  shown 
in  Figure  21  without  ai  being  present. 

B.  Case  (2)  X near  a3  ± an. 

The  governing  equations  for  these  cases  were  developed  in  Sections 
(6.2)  and  (6.3)  and  are 

x o a 

oj  i C i a i a3sinBi  + [ §■  ri222ai  + | 3i2a2]sinY  = 0 

(9.26a) 

(joxC2a2  + siny  = 0 (9.26b) 

01 3 C 3 a 3 + a2sinB2=  0 (9.26c) 

p = 2^rtr C0SY  + felt"  cos8*  (9-26d) 

afcosy  s 2a2  { 7^  ajcosg!  - [!•  ri222a|  + 4 Si2a,.]cosY  } 

c a 4 x (9.26e) 
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where 

Bi 


Y 
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pT2  - 

ar 

- a. 

for 

(JO  3 

+ Wl) 

r 

J 

} (9.27a) 

pT2  + 

ar 

• aJ 

for 

H>3 

- CDll 

- ak 

(9.27b) 

pT2  - 

ar 

- ak 

for 

W3 

+ Wll 

> (9.27c) 

pT2  - 

ar 

+ ak 

for 

UJ3 

- WiJ 

Equations  (9.26)  do  not  readily  admit  a closed-form  solution  so  a 
Newton-Raphson  technique  was  used  to  produce  Figures  22,  23,  and  24  for 
the  case  of  X near  uu  + ui3.  A similar  response  to  that  of  the  last 
section  is  evident.  In  the  presence  of  large  damping,  a!  > a2  and  a3. 
Regions  of  stability  exist  as  well  as  jump  phenomena.  Upon  examining 
the  response,  amplitude  vs  detuning,  Figure  24,  we  note  that  there  is  a 
loss  of  symmetry  from  the  case  of  X near  2uj. 

9.4  Internal  Resonance,  gj3  - 3ioi 

We  now  investigate  the  response  in  the  presence  of  an  internal 
resonance  of  the  form,  003  near  3wi.  For  the  simply  supported  case  this 
is  within  about  10%.  By  a proper  choice  of  the  boundary  conditions,  for 
example  some  form  of  elastic  supports,  it  would  be  possible  to  obtain  a 
relationship  within  much  less  than  10%.  We  now  investigate  the  response 
in  the  presence  of  an  internal  resonance  for  the  case  of  parametric 
resonance: 

(1 ) X near  2<ui 


f 


(2)  no  parametric  resonance 
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A 

Figure  22.  First  modal  amplitude  as  a function  of  aerodynamic 
detuning.  A,  for  constant  frequency  of  excitation, 
p(ep  * X - uj  - W3),  and  constant  amplitude  of 
excitation,  F.  stable, unstable. 
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Figure  24. 
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First  modal  amplitude  as  a function  of  excitation 
frequency,  p(cp  = X - oil  - u)3  ),  for  constant 
amplitude  of  excitation,  F,  and  constant  aerodynamic 
detuning,  A.  stable, unstable. 
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A.  Case  (1 ) Where  X - 2u>i. 

Equations  (9.22)  remain  valid  and  the  entire  discussion,  concerning 
the  response  of  the  first  and  second  modes,  of  Section  9.3  is  valid.  The 
only  change  is  that  the  third  mode  may  now  be  drawn  into  the  response 
through  equations  (7.6c,f)  or 

u)3C3a3  - gi'3  222a2sin8  + -p*2  siny  a 0 (9.28a) 

Gl>  3 d 3 ("2^”  ” O ) + 3 22232C0S8  “ ®2COSY  = 0 (9.28b) 

where  for  e = .001 


= 22.066  - 3(8,1076)  = _ ?1  36 

(.001)  /z 


(9.28c) 


Using  a Newton-Raphson  technique  we  find  that  the  response  is  identical 
to  that  of  Section  9.3  for  no  internal  resonance.  The  relative  sizes 
of  ai  and  a2  are  controlled  by  the  damping  and  detuning  of  the  parametric 
resonance.  The  amplitude  a3  is  very  small  and  the  response  appears  as 
shown  in  Figures  16-21.  To  obtain  a qualitive  feel  for  the  conditions 
under  which  a3  is  drawn  into  the  response,  we  let  the  amplitude  F 
of  the  forcing  term  be  small.  Then  equations  (9.28)  can  be  combined  to 
yield 

r 3222  a| 


a3 


8w3  / C2j  + - o) 


(9.29) 


Thus,  a3  will  be  small  as  long  as  the  radical  of  (9.29)  is  large.  For 
small  damping  the  radical  will  be  small  only  If 
3c 


? 


a ~ 0 


(9.30) 
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Since  the  internal -resonance  detuning,  a,  is  fixed  by  virtue  of  the  form 
of  the  boundary  conditions,  then  the  parametric  detuning  must  be  very 
nearly  2a/ 3 for  a 3 to  be  significant.  However,  in  this  case  a is  large 
and  thus  p must  also  be  large.  But  from  Section  9.3  we  note  that  a 
large  detuning  caused  ai  to  be  so  large  that  a2  was  negligible  by  com- 
parison. So  that  for  large  detuning  of  both  internal  and  parametric 
resonances  the  result  is  essentially  that  of  Figures  16,  17,  and  18. 

As  previously  discussed,  the  detuning  of  the  internal  resonance  is 
governed  by  the  boundary  conditions.  For  small  damping  and  small 
detuning  the  response  is  depicted  essentially  as  described  by  Figures  19 
and  20  with  the  addition  of  the  third  mode;  see  Figures  25  and  26. 

Thus  any  of  the  three  modes  can  be  made  greater  than  the  other  two, 
depending  on  the  type  of  boundary  conditions,  damping,  and  detuning  of 
the  parametric  excitation.  A very  striking  representation  of  the 
relative  sizes  of  the  modes  is  portrayed  in  Figures  27  and  28.  We  note 
here  that  the  response  of  the  third  mode  is  the  same  as  that  for  a 
linear  oscillator  with  a single  degree-of- freedom.  The  second  mode 
amplitude  a2  merely  acts  as  the  external  forcing  function.  A comparison 
with  Figure  18  can  be  made  and  we  see  how  the  decrease  in  damping  allows 
the  second  and  third  modes  to  enter  the  response.  The  second  mode  can 
be  seen  to  be  larger  than  the  first  when  the  detuning  of  the  excitation 
is  near  zero.  Detuning  the  internal  resonance  shifts  the  response  of 
the  third  mode  considerably. 

B.  The  Case  of  No  Parametric  Resonance. 

This  is  essentially  a special  case  of  the  previous  one.  In  the 
absence  of  a parametric  resonance  the  response  for  the  first  two  modes 

1 


Cn*  0.05 
cr  • 0.05 

P • 0.05 
A * 200 


Figure  25.  The  modal  amplitudes  as  functions  of  the  amplitude 

of  excitation,  F,  for  constant  frequency  of  excitation, 

p(ep  * X - 2oji),  constant  aerodynamic  detuninq.  A, 

and  internal  resonance,  a(ca  * <u3  - 3wi). 

stable,  unstable. 


Figure  2 6.  The  modal  amplitudes  as  functions  of  the  aerodynamic 
detuning.  A,  for  constant  frequency  of  excitation, 
a(eo  * X - 2u>i),  constant  amplitude  of  excitation, 

F,  and  internal  resonance  a(ea  = a>j  - 3uu). 

stable.  unstable. 


Figure  28.  The  modal  amplitudes  as  functions  of  the  frequency  of 
excitation,  p(ep  » X - 2<ui),  for  constant  amplitude 
of  excitation,  F,  constant  aerodynamic  detuning,  A and 

internal  resonance  o(ecr  = 013  - 3<m).  stable, 

unstable. 
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remains  the  same,  or 

4aji Ci  C,  ~ A8i  ? 

Vi*  Ti 222 

aj  = 2coiC2a2 

and  the  third  mode  response  is 
1 

a3  = 8 r3.^-a-i (9.35c) 

/ (03  (C3  + 0Z) 

The  response  is  as  shown  in  Figure  29.  Again  any  of  the  three  modes  may 
be  the  largest  depending  on  the  damping  and  detuning  of  the  internal 
resonance.  For  large  damping  &x  > a2  or  a3.  For  small  damping 
a 3 or  a2  >aj  depending  on  the  detuning  of  the  internal  resonance.  In 
addition,  the  larger  the  dynamic  pressure  is  the  more  likely  the  third 
mode  will  be  greatest. 

9.5.  Summary 

The  method  of  solution  described  in  Chapters  5,  6,  and  7 is  used  to 
solve  numerical  examples  of  a simply  supported  plate  undergoing  cylindri 
cal  bending.  Various  combinations  of  internal  and  parametric  resonances 
have  been  examined. 

When  neither  internal  nor  parametric  resonances  are  present  the 
response  depends  on  the  extent  to  which  the  dynamic  pressure  exceeds  the 
critical  value  for  inducing  flutter.  The  amplitude  of  the  response 
increases  as  the  dynamic  pressure  Increases.  Both  modes  associated  with 
the  repeated  frequency  are  excited.  Only  one  stable  solution  exists  and 
it  is  either  trivial  or  nontrivial  depending  on  the  dynamic  pressure. 


(9.35a) 

(9.35b) 
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Figure  29.  The  modal  amplitudes  as  functions  of  the  aerodynamic 


detunings.  A,  for  constant  internal  resonance 
o(co  = uij  - 3wi)  and  no  parametric  resonance. 
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When  a parametric  resonance  of  the  form  X near  2wi  or  id3  ± u>i  Is 
present,  then  two  nontrivial  solutions  can  exist  but  only  one  is  stable. 

The  only  modes  excited  are  associated  with  the  repeated  frequency. 

Jump  phenomona  can  be  produced  by  varying  the  frequency  or  amplitude  of  the 
parametric  excitation  or  the  dynamic  pressure. 

When  both  internal  and  parametric  resonances  are  present  all  three 
modes  can  be  excited.  The  detunings  of  internal  and  parametric  resonances 
have  a great  effect  on  the  third  and  second  modes,  respectively.  When 
the  detuning  of  the  internal  resonance  is  large,  as  in  the  case  of  the 
simply  supported  plate  considered,  we  see  that  the  response  is  essentially 
identical  to  the  case  of  no  internal  resonance.  Jump  phenomena  exist 
regardless  of  which  amplitude  is  largest. 
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APPENDIX  A 


SIMILARITY  TRANSFORMATIONS 

A. 1 General 

Given  a nonsingular  matrix  Z,  then 

[J]  - [Z_1][B][Z]  (Al) 

is  a similarity  transformation  where  the  matrices  [J]  and  [B]  have  the 
same  eigenvalues.  The  choice  of  [Z]  can  determine  the  complexity  of  the 
form  of  [J],  the  simplist  being  the  diagonal  form  when  all  the  eigen- 
values as  distinct.  If  the  eigenvalues  are  not  all  distinct  [J]  may 
still  be  diagonal,  if  the  eigenvectors  corresponding  to  the  repeated 
roots  are  distinct.  When  the  eigenvectors  corresponding  to  repeated 
roots  are  identical,  then  the  matrix  can  at  best  take  on  the  Jordan 
form,  where  the  eigenvalues  appear  as  diagonal  terms  and  a 1 (one) 
appears  as  an  off-diagonal  term  adjacent  to  the  repeated  roots;  see 
equation  (9.11b).  The  number  of  off-diagonal  terms  depends  on  the 
number  of  repeated  frequencies  or  more  precisely  the  number  of  times 
a given  root  is  repeated,  and  the  corresponding  number  of  independent 
eigenvectors.  The  latter  depends  on  the  rank  of  the  matrix  which  is  not 
the  major  topic  of  discussion  of  this  appendix.  However,  a brief  aside 
is  In  order.  The  reader  must  accept  what  follows  In  this  example  with 
faith  or  confirm  if  with  rigor  in  a mathematics  text  such  as  Kreyszig 
(1972)  (pp.  256-279). 

Given  an  n x n matrix  with  one  repeated  frequency,  then.  If  the 
rank  of  the  matrix  is  r * n - 1,  there  is  only  one  linearly  Independent 
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eigenvector  associated  with  the  repeated  frequency.  If  the  rank  is 
r = n - 2 then  there  are  two  linearly  independent  elgnevectors  associated 
with  the  repeated  frequency.  The  last  case  will  allow  generation  of  a 
diagonal  matrix  while  the  first  will  produce  the  Jordan  form  with  one 
off  diagonal  term. 

The  problem  is  one  of  determining  the  matrix  [Z]  which  will  put  [B] 
in  Jordan  form.  Mathematics  texts  indicate  [Z]  exists  and  give  examples 
to  prove  it  or  direct  the  student  to  verify  a sample  problem  with 
mechanical  computations.  However,  nowhere,  does  there  appear  to  be  a 
method  to  determine  the  matrix  [Z]  in  the  first  place.  That  Is  precisely 
the  purpose  of  this  Appendix. 

Finding  the  eigenvalues  and  eigenvectors  of  [B]  is  a relatively 
trivial  problem  with  the  advent  of  sophisticated  library  Computer 
programs.  Using  that  as  a starting  point  we  proceed  with  a Newton- 
Raphson  iterative  technique  to  generate  the  matrix  [Z]  from  the  set  of 
eigenvectors. 

A. 2 Determining  the  Form  of  [J], 

We  start  with  the  matrix  [B]  and  use  any  standard  library  computer 
program  which  will  return  eigenvalues  and  eigenvectors.  An  Inspection 
of  eigenvalues  will  Indicate  immediately  if  any  are  repeated.  If  they 
are  too  numerous,  a simple  sort  routine  can  be  employed  to  order  them 
from  highest  to  lowest  to  aid  in  Identifying  repeated  roots.  Once  the 
repeated  roots  are  identified,  the  corresponding  eigenvectors  are 
compared.  The  best  way  Is  to  normalize  them  by  dividing  the  elements, 
usually  returned  as  a column  vector,  by  the  value  of  the  first  element. 
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The  vectors  which  remain  are  either  identical  or  they  are  not.  If  they 
are  identical  then  the  Jordan  form  will  evolve.  There  need  not  be 
any  discussion  of  rank,  linear  independence  or  whatever. 

A. 3 Newton-Raphson  Iteration  Technique 

The  Newton-Raphson  technique  is  especially  useful  in  solving 
nonlinear  problems.  It  is  "initial  guess"  sensitive  to  some  degree  but 
converges  rapidly  if  the  system  is  fairly  well  behaved  near  the  solution. 
The  basic  problem  is  one  where  we  are  required  to  find  the  values  of  a 
vector  {x}  which  will  make 

{f ( x ) } = 0 (A2) 

For  a 2 x 2 example  we  expand  in  a Taylor  series  and  obtain 

fi(xi)  = fi(xio)  + (xio)Axi  + (xio)Ax2  + ...  (A3a) 

f 2 (x2 ) = f 2 (x2 o ) + 1^-  (x2„)Axi  + 1^-  (x20)Ax2  + ...  (A3b) 

where  the  x„„  are  the  initial  guesses.  The  Ax  are  the  increments 
which  will  cause  all  the  f i (xp ) to  be  zero  simultaneously,  so  we  write 

fxi"  AXl  + "fx^"  AXz  = " (A4a ) 

Ixf  AXl  + lx?"  A*2  * ” f2(^2o)  (A4b) 

and  we  have  a linear  system  which  can  be  solved  for  the  Axn  to  add  to 
the  xn(J  and  proceed  to  another  iteration  until  we  reach  some  pre- 
determined limit  of  convergence  in  either  the  f^  or  Ax^ . 
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A.  4 The  [2]  Matrix  With  One  Repeated  Frequency 


The  [Z]  matrix  is  constructed  as  follows.  The  eigenvectors  are 
sorted  so  that  the  repeated  vectors  are  in  the  first  two  columns  and 
are  normalized  so  the  top  element  in  each  column  is  one.  We  seek  a 
column  vector  to  replace  the  first  column,  the  repeated  eigenvector. 
The  vector  function,  f,  represents  the  first  column  of  the  matrix  [J] 
below  the  first  element;  that  is 


{Ji>  = 


1 

Mxi) 

fz(x2) 


Fn-l(xnl'; 


(A5a) 


We  wish  the  final  result  to  be 


(A5b) 


L°J 

If  the  functional  form  f(xn)  were  known,  we  would  have  no  problem  in 
evaluating  the  partial  derivatives 


3f 

3x 


(A6) 
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Here  we  must  construct  them  numerically.  We  begin  by  assuming  the  first 
column  of  the  [Z]  matrix  to  be 


Then  we  compute 

[J]  = [Z_1][B][Z]  (A8) 

The  first  column  of  [J]  becomes 
1 

i(x0) 

2 ( X o ) 
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and  so  on,  until  we  have  a system  of  the  form  of  (A4)  and  we  write  for 
r = n - 1 


where  the  right  side  was  computed  from 


and  the  factor  of  one  applied  to  the  first  element  to  satisfy  (A5b). 
Equations  (A15)  may  now  be  solved  for  the  Axi  and  added  to  one  to  form  a 
new  guess  for  the  x^ . Now  the  product  of  (Al)  is  computed  tc  determine 
if  the  Jordan  form  has  been  achieved.  Convergence  Is  achieved  by 
specifying  the  accuracy  with  which  (A17)  is  satisfied: 

[Z_1][B][Z]  - [J]  : 0 (A17) 

In  this  problem  one  step  was  required  to  produce  an  accuracy  of  10" 10  in 
the  element  using  double  precision. 
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PARAMETRICALLY  EXCITED  NONLINEAR 


MULTI-DEGREE-OF-FREEDOM  SYSTEMS 
by 

Edward  G.  Tezak 
(ABSTRACT) 

An  analysis  of  parametrically  excited  nonlinear  multi -degree-of- 
freedom  systems  is  presented.  The  nonlinearity  considered  is  cubic  and 
small  so  that  the  system  of  equations  is  weakly  nonlinear.  Modal  damping 
is  included  and  the  parametric  excitation  is  harmonic.  The  systems 
examined  include  those  with  distinct  natural  frequencies  as  well  as 
those  with  a single  repeated  frequency.  The  significant  role  played  by 
the  existence  of  an  internal  resonance  is  explored  in  depth.  The  de- 
rivative-expansion version  of  the  method  of  multiple  scales,  a per- 
turbation technique,  is  used  to  develop  solvability  conditions  for  the 
various  combinations  of  internal  and  parametric  resonances  considered. 

Regions  where  trivial  and  nontrivial  solutions  exist  are  defined 
and  the  stability  of  the  solutions  within  each  region  is  discussed. 
Nontrivial,  unstable  solutions  have  been  shown  to  exist  in  regions  where 
nontrivial  stable  solutions  are  known.  Numerical  solutions  do  not  hint 
at  the  existence  of  these  solutions. 

The  role  of  internal  resonance  in  parametrically  excited  systems  is 
explored.  Strong  modal  interaction  is  demonstrated  as  a consequence  of 
the  presence  of  the  cubic  nonlinearity  and  the  Internal  resonance. 
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Because  of  this  modal  coupling,  modes  other  than  the  one  excited  can 
dominate  the  response.  A multiplicity  of  jumps  is  shown  to  exist. 


Parametric  excitation  in  the  nonlinear  flutter  problem  has  been 
examined  in  detail  for  the  first  time.  The  effect  of  the  parametric 
excitation  can  raise  the  flutter  speed  if  the  amplitude  is  small  and 
lower  it  if  the  amplitude  is  large. 

Limit-cycle  behavior  in  the  flutter  problem  is  developed  in  a 
unique  analytic  way.  The  condition  which  predicts  the  onset  of  flutter 
is  developed  by  using  a linear  analysis.  The  solution  grows  without 
bound.  Interestingly,  the  same  condition  in  the  nonlinear  analysis 
predicts  the  existence  of  a real  nontrivial  solution  with  a finite 
amplitude,  the  so-called  limit-cycle. 


